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Abstract. A third workshop on small- a; physics, within the Small-x Collaboration, was held in Hamburg 
in May 2004 with the aim of overviewing recent theoretical progress in this area and summarizing the 
experimental status. 
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very encouraging although still very far from the precision 
in QED. 

However, there is a domain, stiU in the perturbative 
regime, where our understanding is lacking. This is the 
region of high energy and moderate momentum transfer, 
such as small-x Deeply Inelastic Scattering (DIS) as mea- 
sured at HERA and low to medium E± jet production 
at the Tevatron. In this region, the collinear factorization 
must break down as the perturbative expansion becomes 
plagued by large logarithms of the ratio between the total 
collision energy and the momentum transfer of the hard 
sub-process, which needs to be resummed to all orders to 
obtain precision predictions from QCD. These logarithms 
arise from the large increase of the phase space available 
for additional gluon emissions, resulting in a rapid rise 
of the gluon density in hadrons with increasing collision 
energy or, equivalently, decreasing momentum fraction, x. 

In this high energy limit, QCD is believed to be cor- 
rectly approximated by the BFKL evolution [6-8], and 
cross sections should be possible to predict using k±- 
factorization [9-12] where off-shell matrix elements are 
convoluted with unintegrated parton densities obeying 
BFKL evolution. However, so far the precision in the pre- 
dictions from fcj^-factorization has been very poor. Al- 
though BFKL evolution correctly predicted the strong rise 
of the F2 structure function with decreasing x at HERA on 
a qualitative level, it turned out that the next-to-leading 
order corrections to BFKL are huge [13,14], basically mak- 
ing any calculation with leading-logarithmic accuracy in 
fc_L-factorization useless. 

Several attempts have been made to tame the NLO 
corrections to BFKL by e.g. matching to the collinear 
limit [15] and matching this with off-shell matrix elements 
or impact factors calculated to NLO. Another strategy is 
based on the fact that a large part of the NLO corrections 
to BFKL can be traced to the lack of energy and mo- 
mentum conservation in the LO evolution [16]. Although 
energy and momentum is still not conserved in NLO 
evolution, the contributions from ladders which violates 
energy-momentum conservations are reduced. Amending 
the leading-logarithmic evolution with kinematical con- 
straints, either approximately in analytical calculations 
[17] or exactly in Monte-Carlo programs [18-21], should 
possibly lead to more reasonable QCD predictions, al- 
though still formally only to leading logarithmic accuracy. 
However, so far none of these strategies have been able to 
fulfill their ambitions, and the reproduction of available 
data is still not satisfactory. 

The plot thickens further when considering the in- 
crease in gluon density at small x. At high enough en- 
ergy the density of gluons becomes so high that they must 
start to overlap and recombine, and we will encounter the 
phenomena of multiple interactions, saturation and rapid- 
ity gaps. In the non-perturbative region these phenomena 
have already been established, but there is currently no 
consensus on whether effects of recombination of pertur- 
bative gluons have been seen at e.g. HERA. Perturbative 
recombination would require non-linear evolution equa- 
tions, which then also could break fcj^-factorization. 



In our first review [22] we focused on the theoretical 
and phenomenological aspects of fc^ -factorization, while in 
the second [23] we also gave an overview of experimental 
results in the small- a; region. In this third review we will 
continue to present recent developments in these areas, 
but also give an overview and introduction to saturation 
effects and non-linear evolution. 

The layout of this report is as follows. First we dis- 
cuss some recent developments of ^^-factorization in sec- 
tion 13 starting with the unintegrated parton densities 
(section 12. 2|) and doubly unintegrated parton densities 
(|2.3(l and continuing with recent advances in NLO cal- 
culations 1)2.41 and I2.6|l . Then, in section 13 we describe 
some phenomenological applications of fcj^-factorization, 
looking at how to use them to obtain QCD predictions for 
heavy quark (|3.1I) and quarkonium H3.4|l production. In 
section^ we present the recent investigations by March- 
esini and Mueller relating some aspects of jet physics to 
BFKL dynamics, which could make it possible to study 
this kind of evolution also in other environments. In sec- 
tion [3 we give an introduction and overview of saturation 
phenomena and non-linear evolution. Section |S| also deals 
with saturation, but in the context of the so-called AGK 
cutting rules which enables us to relate saturation with 
multiple scatterings and diffraction. In section|7|we review 
some recent experimental results relating to the issues in 
the previous sections, beginning with multiple interactions 
and underlying events in section [V. II followed by rapidity 
gaps between jets in 17.21 jet-production at small- a; in 17.31 
and production of strange particles in DIS in section 17.41 
Finally we present a brief summary and outlook in sec- 
tion |H1 



2 The fcx-factorization formalism 

Main author H. Jung 

In the high energy limit, cross sections can be cal- 
culated using k± -factorization [9-12] with convolution 
of a off-shell {k± dependent) partonic cross section 
(T(-,fc^) and an k± - unintegrated parton density func- 
tion :P(z, A:^): 



a= j—dH.^a{-,kl)T{z,k\ 

J Z Z 



(1) 



The unintegrated gluon density J^(z,fc^) is described by 
the BFKL [6-8] evolution equation in the region of asymp- 
totically large energies (small x). An appropriate descrip- 
tion valid for both small and large x is given by the CCFM 
evolution equation [24-27], resulting in an unintegrated 
gluon density, Aix, fc^, g^), which is a function also of the 
additional scale, q. Here and in the following we use the fol- 
lowing classification scheme: xQix^k^) describes DGLAP 
type unintegrated gluon distributions, xT(x,k\^ is used 
for pure BFKL and xA(x, k^,q^) stands for a CCFM type 
or any other type having two scales involved. Different 
approaches to the unintegrated parton density functions 
have been discussed in detail in [22,23]. 
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While still being formally at leading order, the uninte- 
grated gluon densities incorporate effects from the next- 
to-leading order in the coUinear approach [28] . This is dis- 
cussed in more detail in the next subsections. To further 
connect to the uncertainty estimates of cross section cal- 
culated in the coUinear approach, the change of the renor- 
malization and factorization scales are used to estimate 
the influence and size of higher order corrections. In [29] 
the CCFM unintegrated PDFs are determined such that 
the structure function F2 as measured at HI [30,31] and 
ZEUS [32,33] can be described after convolution with the 
off-shell matrix element. This fit is repeated for the renor- 
malization scale in the off-shell matrix element varied by a 
factor of 2 up and down, resulting in new sets of PDFs [29] , 
set A0+ and set A0-. These PDFs are compared with the 
central set set AO in Fig.Q] 



charm anti-charm pair. In [44], also an alternative to this 
was discussed, namely to measure the quantity 



/bmax>^Pmi„;W^) 



(t{W) 



(3) 



where p^a^^ = Taa:K{dpl^,dpl^), p^i„ = min(dp? 4, dpiJ 
and /c is a constant. This quantity would be a measure of 
the spread in the p\^ x p\ ^ plane. Yet another possibility 
would be a direct reconstruction of Xg and k^ ^ from (DIS) 
multijet events, thereby mapping the unintegrated gluon 
density directly. 

The unintegrated gluon density could also be con- 
strained from global fits. So far, only fits to F2 have been 
made [45], and a global fit using various data such as 
forward jets, 2-|-n jets, heavy quarks and azimuthal jet- 
jet correlations would further constrain the unintegrated 
gluon density. 




Fig. 1. Comparison of the CCFM uPDF obtained after chang- 
ing the renormalization scale in the off-shell matrix element by 
a factor 2 up and down. 



2.1 Future fits of uPDF parameterizations 

Main author M. Hansson 

There are a number of possible measurements sensi- 
tive to the transverse momentum of the propagating glu- 
ons in the gluon ladder, and thereby suitable for investi- 
gations concerning the unintegrated gluon density of the 
proton. One possible observable is the difference in az- 
imuthal angle, A(t)*, of a dijet system in the hadronic cen- 



ter of mass frame. The differential cross section 



dA<p* 



has 



been measured at the Tevatron [34-39] and only recently 
at HERA [40,41]. The quantity 



S^ 



j;^Nd,,et{Ar,x,Q^)dAr 

^ Nd^Jet{M\X:Q'')dA4>-' 



(2) 



first proposed in [42] , has been measured [43] and showed 
a large sensitivity to the unintegrated gluon density. An- 
other measurement, proposed in [44], would be to mea- 
sure , a'^j 2 where dpj , are the transverse momenta of a 



2.2 The need for doubly unintegrated parton density 
functions 

Main author J. Collins 

Conventional parton densities are defined in terms of 
an integral over all transverse momentum and virtuality 
for a parton that initiates a hard scattering. While such 
a definition of an integrated parton density is appropriate 
for very inclusive quantities, such as the ordinary struc- 
ture functions Fi and F2 in DIS, the definition becomes 
increasingly unsuitable as one studies less inclusive cross 
sections. Associated with the use of integrated parton den- 
sities are approximations on parton kinematics that can 
readily lead to unphysical cross sections when enough de- 
tails of the final state are investigated. 

We propose that it is important to the future use of 
pQCD that a systematic program be undertaken to re- 
formulate factorization results in terms of fully uninte- 
grated densities, which are differential in both transverse 
momentum and virtuality. These densities are called "dou- 
bly unintegrated parton densities" by Watt, Martin and 
Ryskin [46,47] (discussed in the next section), and "par- 
ton correlation functions" by Collins and Zu [48]; these 
authors have presented the reasoning for the inadequacy, 
in different contexts, of the more conventional approach. 
The new methods have their motivation in contexts such 
as Monte-Carlo event generators where final-state kine- 
matics are studied in detail. Even so, a systematic refor- 
mulation for other processes to use unintegrated densities 
would present a unified methodology. 

These methods form an extension of fcj^-factorization, 
which has so far been applied in small- a; processes and, 
as the CSS formalism [49], in the transverse-momentum 
distribution of the Drell-Yan and related processes. 

The problem that is addressed is nicely illustrated by 
considering photoproduction of cc pairs. In Figs. |21 we 
compare three methods of calculation carried out within 
the Cascade event generator [20,50]: 
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Fig. 2. (a) and (b): Comparison between use of simple LO 
parton model approximation and of the use of k± densities for 
the p± of cc pairs in photoproduction, and for the x^. (c) and 
(d): Comparison of use of k± densities and full simulation. 



shown in plot (b), x-y stays close to unity. Neglecting the 
full recoil mass m is equivalent of taking k — tj— ^ with 
k^ being the virtuality of the gluon, fc^ its transverse mo- 
mentum and X its light cone energy fraction. This gives 
a particular value to the gluon's k~ . When we also take 
into account the correct virtuality of the gluon, there is no 
noticeable change in the pxce distribution — see Fig.|2Ic) 
(dotted line) — since that is already made broad by the 
transverse momentum of the gluon. But the gluon's k~ is 
able to spread out the x-y distribution, as in Fig.[5fd) with 
the dotted line. This is equivalent with a proper treatment 

of the kinematics and results in fc^ = — -^_^™ , which 
can be significant for finite x. Clearly, the use of the sim- 
ple parton-model kinematic approximation gives unphys- 
ically narrow distributions. The correct physical situation 
is that the gluon surely has a distribution in transverse 
momentum and virtuality, and for the considered cross 
sections neglect of parton transverse momentum and vir- 
tuality leads to wrong results. It is clearly better to have 
a correct starting point even at LO, for differential cross 
sections such as we have plotted. 

Therefore it is highly desirable to reformulate pertur- 
bative QCD methods in terms of doubly unintegrated par- 
ton densities from the beginning. A full implementation 
will be able to use the full power of calculations at NLO 
and beyond. 



— Use of a conventional gluon density that is a function 
of parton x alone. 

— Use of a k± density that is a function of parton x and 
k± . These are the objects usually called "unintegrated 
parton densities" . 

— Use of a "doubly unintegrated density" that is a func- 
tion of parton a;, k± and virtuality, that is, of the com- 
plete parton 4-momentum, in Cascade taken after the 
full simulation of the initial state parton showering. 

The partonic subprocess in all cases is the lowest order 
photon-gluon-fusion process j + g — > c -\- c. Two differ- 
ential cross sections are plotted: one as a function of the 
transverse momentum of the cc pair, and the other as a 
function of the x-y of the pair. By x-y is meant the frac- 
tional momentum of the photon carried by the cc pair, 
calculated in the light-front sense as 

Here i^e is the electron beam energy and the coordinates 
are oriented so that the electron and proton beams are in 
the —z and -\-z directions respectively. 

In the normal parton model approximation for the 
hard scattering, the gluon is assigned zero transverse mo- 
mentum and virtuality, so that the cross section is re- 
stricted to "Ptcc = and x-y = 1, as shown by the solid 
lines in Fig. |5fa,b). When a k± dependent gluon density 
is used, quite large gluonic k± can be generated, so that 
the ptcc distribution is spread out in a much more phys- 
ical way, as given by the dashed line in Fig.[2Ia). But as 



2.3 Doubly unintegrated PDFs 

Main author G. Watt 

The notation for the two-scale unintegrated gluon dis- 
tribution, xA{x,k'^,q^), used in [22,23] and elsewhere in 
this report, is related to that used in this section by 

xA{x,k'i,f) ^ fg{x,kt,fi'^)/kt. 



(4) 



2.3.1 Unintegrated PDFs from integrated ones 



Existing analyses of the CCFM equation are based on 
numerical solution via Monte Carlo methods. Kimber, 
Martin and Ryskin [51] showed that, in a certain ap- 
proximation, it is possible to obtain two-scale UPDFs, 
fa{x,k1,y?), from single-scale distributions, with the de- 
pendence on the second scale /i introduced only in the last 
step of the evolution. It was found that this "last-step" 
prescription gave similar results whether the single-scale 
distributions were evolved with a unified BFKL-DGLAP 
equation [52] or purely with the DGLAP equation, indicat- 
ing that angular ordering is more important than small- 
X effects. Here, we summarize the procedure [46, 51] for 
obtaining UPDFs from the conventional DGLAP-evolved 
integrated PDFs, a{x,fj?) = xg{x,fj?) or xq{x,^'^). 

The UPDFs are constructed to satisfy the normaliza- 
tion conditions 



'^t 



a{x,f?), 



(5) 
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which are ensured by defining the UPDFs to be [46,51] 
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Hx.z,k^^.^^) 
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b=g,q 



dzP, 



>a,{z)b[-^,k^), (6) 



where the Sudakov form factors are 



Fig. 3. Illustration of (z, fct)-factorization for the doubly- 
unintegrated gluon distribution, fg{x,z,kt , fi^), shown in the 
final diagram. In the first two diagrams the penultimate par- 
ton in the DGLAP evolution chain, with 4-momentum fc„_i — 
{x/z)p, splits into a gluon with 4-momentum kn = k — 
xp~ I3q' + k±. 



exp 



'dK? asinj) 



fe? «;? 



27r 



b=g,q 



dC CAa(C), (7) 2.3.2 Doubly-unintegrated PDFs 



and Pba are the unregulated LO DGLAP splitting kernels. 

In addition, it is necessary to apply angular-ordering 
constraints due to color coherence, which regulate the sin- 
gularities in © and Q arising from soft gluon emission. 
These constraints are not applied for quark emission where 
there is no "coherence" effect. The explicit expressions for 
the unintegrated gluon and quark distributions are given 
in [46]. 

This approach to UPDFs amounts to relaxing the 
DGLAP approximation of strongly-ordered transverse 
momenta along the evolution chain only in the last evo- 
lution step. If we consider DIS in the Breit frame, where 
the proton has 4-momentum p and the virtual photon has 
4-momentum q, then the penultimate parton in the evo- 
lution chain, with 4-momentum fc„_i = (x/z)p, splits to 
a final parton with 4-momentum 



kn = k — {k^ , k , kt) = X p — f3 q' + k±, 



(8) 



where the plus and minus components are k = A:° ± fc"^. 
In the Breit frame: 



P= {Q/xB], 0,0), 

q' = q + XBjP= (0,(3,0), 

fci = (0,0,fct), 



(9) 
(10) 
(11) 



/2 2 

q , q 



-Q^ and fc^ 



fc|. The 



so that p"^ = 

condition that the parton emitted in the last evolution 

step is on-shell, {k„-i — fcn)^ = 0, gives 



P 



XBj 



kl 



X (l-z)Q2' 



(12) 



so fc2 = — /c^/(l — z). In the high-energy (small-x) limit, 
where gluons dominate, we have z — > 0, so A: ~ a;p -I- fc^ 
and k^ ~ —kf- Cross sections can then be calculated using 
the /ct-factorization formalism. 



1 dx 

X JQ 



1.2 Jgy 



— I ^h{x,klt?)a^^ , (13) 

Bj X Jo '^t 



where the partonic cross section a^ ^ is calculated with 
an off-shell incoming gluon. 



Away from the high-energy limit, where we have finite 
z, the partonic cross section of (|13|l will necessarily have 
some z dependence through the q' component, i.e. the 
minus component, of the 4-momentum k (jHll. Therefore, 
we should consider doufeZj/- unintegrated PDFs (DUPDFs), 
fa{x, z,kf, fM^), which satisfy 

1 

dzfa{x,z,kl,^i^) = fa{x,k'^,p,'^). (14) 

From ®, the DUPDFs are 



fa{x,Z,kt,^i ) ^Ta{h,fl ) 



2tt 



J2Pabiz)b[^,k^), (15) 



b=g,q 

apart from the angular-ordering constraints. The explicit 
expressions for the doubly-unintegrated gluon and quark 
distributions are given in [46]. The /c^ -factorization for- 
mula (|13|l is then generalized to the "(z, A:i)-factorization" 
formula [46] 



a'' P = 

Y. f - fdz r^/,(.,.,fc,^^2)^.-a*^ (ig) 
a=g,q''^Bj ^ Jx Jo «t 

Note that fa{x,z,k^,iJ?') are linear densities in z, but log- 
arithmic in X and fc^. This idea is illustrated in Fig.|21for 
the case a = g. It is not immediately obvious how the 
partonic cross sections ct^ " in (|16|l should be calculated. 
Recall that they can be written 



a= fd^\M\'^/F, 



(17) 



where d^ is the phase space element, [A^l^ is the squared 
matrix element, and F is the flux factor. The phase space 
element d'P can be calculated with the full kinematics, 
that is, with k = xp — f3q' + k±. The flux factor F is 
taken to be the same as in collinear factorization (and 
in /ct-factorization) , that is, F ~ Axp ■ q. The last evo- 
lution steps in Fig. 13 only factorize from the rest of the 
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Fig. 4. Comparison with HI inclusive jet production data [53] 
in three pseudorapidity (r;'"*°) bins. The predictions of the 
(2, fct)-factorization approach based on DUPDFs [46] (which 
is much simpler to implement) are in good agreement with the 
conventional QCD approach. In some bins the predictions of 
the latter approach are hidden beneath the bold lines of the 
(2, fct)-factorization approach, at the respective order. 



diagram, to give the LO DGLAP splitting kernels, in the 
leading logarithmic approximation (LLA), that is, in ei- 
ther the coUinear (kt -^ 0) or high-energy [z -^ 0) limits. 
Therefore, \M}\' should be evaluated with either k = xp 
or fc = a; p + fc^ , in order to provide the factorization be- 
tween the DUPDF and the subprocess labeled a in Fig. 13 
For the specific case of inclusive jet production in DIS and 
working in an axial gluon gauge, it was observed in [46] 
that the main effect of the "beyond LLA" terms (pro- 
portional to (3 dJ) was to suppress soft gluon emission, 
and that these terms made a negligible difference to the 
cross section when the angular-ordering constraints were 
applied. 

The prescription adopted in [46] was to evaluate |A^|^ 
in the coUinear approximation (fc = xp), so that a (z, fcj)- 
factorization calculation approximately reproduces the 
coUinear factorization calculation starting one rung down 
as in the first two diagrams of Fig. |2| that is, where the 
subprocess is evaluated at one order higher in 0:5 . This was 
demonstrated for inclusive jet production in DIS, where 
the LO subprocess is simply 7*9* — > q. Similarly, a "NLO" 
calculation, where the subprocesses are 7*5* — > qq and 
7*g* -^ qg^ was found to give results close to the conven- 
tional NLO QCD calculation, where the subprocesses are 
0(a|); see Fig. H 

In [47], the (z, A;t)-factorization formalism was ex- 
tended to hadron-hadron collisions and applied to pre- 
dict the pt distributions of vector bosons {V — W, Z) and 
Standard Model Higgs bosons (H). For pT ^ My.Hi fixed- 
order coUinear factorization calculations diverge, with 
hi{Mv,H / Pt) terms appearing in the perturbation series 
due to soft and coUinear gluon emission. Traditional calcu- 
lations combine fixed-order perturbation theory at high px 
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Fig. 5. Pt distribution of W bosons produced at the Teva- 
tron calculated using [z, fct)-factorization [47], compared to D0 
data [56]. 



with either analytic resummation or numerical DGLAP- 
based parton shower formalisms at low pr, with some 
matching criterion to decide when to switch between the 
two. It has been shown in [54, 55] that UPDFs obtained 
from an approximate solution of the CCFM evolution 
equation embody the conventional soft gluon resumma- 
tion formulae. In the framework of (z, A:t)-factorization, 
the lowest order subprocesses are simply qlq2 ^> V and 
9i92 ^ H. A good description was obtained in [47] of 
the Pt distributions of W and Z bosons produced at the 
Tevatron Run 1 over the whole pr range; see Fig. El The 
predicted Higgs pr distribution at the LHC was found to 
reproduce, to a fair degree, the predictions of more elabo- 
rate theoretical studies [57] , in particular the NNLL-I-NLO 
resummation approach of Grazzini et al. [58]; see Fig. 1^1 
Alternative predictions for Higgs production at the LHC 
using the fcj-factorization approach have been made in 
[54,59-61]. 

Note that matrix-element corrections are necessary 
in DGLAP-based parton shower simulations at large px- 
Without such corrections, the herwig parton shower pre- 
diction falls off dramatically at large pr ^ Mh [62]; see 
Fig. El The same effect is observed in herwig predic- 
tions for the Pt distributions of W and Z bosons [63], 
whereas in Fig. the Tevatron data at large pT > Mw 
are well-described without explicit matrix-element correc- 
tions. Also, the (z, fcf)-factorization prediction for Higgs 
production is found to be close to the NLO fixed-order 
result at large pt, see Fig. El suggesting that a large part 
of the subleading terms are included by accounting for the 
precise kinematics in the gl g2 —^ H subprocess. 

The integrated PDFs used as input in [46,47] were de- 
termined from a global fit to data using the conventional 
coUinear approximation [64]. A more precise treatment 
would determine the integrated PDFs, used as input to 
the last evolution step, from a new global fit to data using 
the (z, fcf)-factorization formalism. 
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Fig. 6. pt distribution of SM Higgs bosons produced at the 
LHC with mass 125 GeV calculated using (z, fci)-factorization 
[47] , compared to various resummed and parton shower predic- 
tions which are all matched to fixed-order calculations at large 
Pt (apart from herwig) [57]. 



2.4 NLO BFKL 



Main author J. Andersen and A. Sabio-Vera 



Since the completion of the calculation of the next-to- 
leading (NLL) corrections to the BFKL equation [13, 14] 
for the forward kernel there has been a large activity fo- 
cused on the study of the fundamental properties of the 
NLL gluon Green's function in the Regge limit of QCD at 
high energies [15,65-82]. Recently, a powerful approach 
has been developed which allows for the complete and ex- 
act analysis of the solution at NLL. In Ref. [83] it was 
demonstrated how it is possible to use D = 4 -I- 2e di- 
mensional regularization together with an effective gluon 
mass (A) to explicitly show the cancellation of simple and 
double poles in e. This procedure carries a logarithmic 
dependence in A which numerically cancels out when the 
full NLL BFKL evolution is taken into account for a given 
center-of-mass energy, this being a natural consequence 
of the infrared finiteness of the full kernel. The basis of 
this approach is the iterated form of the solution for the 
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where the strong ordering in longitudinal components of 
the parton emission is encoded in the nested integrals in 
rapidity with an upper limit set by the logarithm of the 
total energy in the process, j/o = Y. The Reggeized form 
of the gluon propagators in the i-channel, uj^ (q), in this 
approach reads 
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with 



^^ ' 4 V3 3 3iVc iVc m2 



(20) 



being the corresponding part in the real emission kernel. 
To complete the real part of the NLL kernel there are other 
more complicated terms in K-r which do not generate e 
singularities when integrated over the full phase space of 
the emissions, for details see Ref. [83]. 

The numerical implementation and analysis of the 
form of solution as in Eq. H18|) was carried out in Ref. [84] . 
At the light of this study the known feature of a lower in- 
tercept at NLL with respect to leading-order (LL) was 
confirmed. As in this approach it is not needed to expand 
on any eigenfunctions there are no instabilities in the en- 
ergy growth. This is highlighted at the left hand side of 
Fig. where the bands correspond to uncertainties in the 
choice of renormalization scale. 

However, the space where the convergence of the per- 
turbative expansion is poor is not in energy but in trans- 
verse momenta. In particular, when the two transverse 
scales entering the forward gluon Green's function are 
of comparable magnitude then the NLL corrections are 
smaller when compared to LL, this can be seen in the bot- 
tom plot of Fig [3 However when the ratio between these 
scales largely departs from unity then the [NLL — LL[ dif- 
ference becomes large, driving, as it is well-known, the 
gluon Green's function into an oscillatory behavior with 
negative values. 
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Fig. 7. Analysis of the gluon Green's function as obtained Fig. 8. Distribution in the number of iterations and angular 
from the NLL BFKL equation. dependence of the NLL gluon Green's function. 



The main advantage of the method here described 
is that the Green's function is generated integrating the 
phase space using a Monte Carlo samphng of the different 
parton configurations. This feature allows for a full con- 
trol of the average multiplicities and angular dependences. 
The former can be extracted from the Poisson-like distri- 
bution in the number of rungs, or iterations of the kernel, 
needed to reach a convergent solution. This is obtained 
numerically in the upper part of Fig. |S1 where we see e.g. 
that for y = 5 it is should be enough to include ^ 15 
rungs/iterations. At the lower part of the same figure the 
angular correlations in the azimuthal angle of dijets with 
similar and large transverse energy, and low hadronic ac- 
tivity in between, is studied in a toy cross-section with 
simplified impact factors. The increase of the angular cor- 
relation when the NLL terms are included in such observ- 
able is a characteristic feature of these corrections. This 
study is possible within this approach in an immediate 
manner because the NLL kernel is treated in full, with- 
out angular averaging, so there is no need to use a Fourier 
expansion in angular variables via the introduction of con- 
formal spins. 



An interesting theoretical development in the context 
of NLL BFKL was the calculation of the forward NLL ker- 
nel in the conformally invariant N = 4 super Yang-Mills 
theory [85,86]. In such field theory the coupling remains 
a constant even at NLL, opening the possibility of finding 
the solution of the BFKL equation in a straightforward 
way because the LL eigenfunctions are also so at NLL. 
In particular, the kernel was calculated for all conformal 
spins in Ref. [85,86] allowing for the direct test of the angu- 
lar structure of the solution as obtained from the method 
here described. This comparison between both approaches 
was performed in Ref. [87]. In this case the gluon Regge 
trajectory reads (with a denoting the coupling constant) 



1 ^\«' 



y-^-?.-c(3) 



(21) 



'^o (q) 



and ^ = a + 

mic dependence. For a precise determination of the contri- 
bution to the gluon Green's function stemming from the 



-j^ — Y2 ) is a constant without logarith- 
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different Fourier components in the aziniuthal angle, i.e. 
/(ka,kb,Y)= J2 ^(|ka|,|kb|,Y)e^"^ (22) 

n— — oo 

it is enough to extract the coefficients of the expansion, 
either using the kernel calculated in [85,86] 



/„(|ka|,|kb|,Y) 



1 



Trlkallkbl 



d-l ^k|y %^„(a,7)Y 



(23) 



27ri Vk^ 



or making use of the iterative solution explained in this 
section [87]: 



fn (jka 



|kb|,Y) = 



y''"|^/(ka,kb,Y)cOs( 



n9). (24) 



The results from these two independent alternatives are 
shown to coincide in Fig. IHl In the upper part the n = 
Fourier component clearly dominates at large energies, de- 
creasing the angular correlations as the energy increases. 
In the lower part it is shown how the convergence in the 
angular variable on the transverse plane is achieved after 
only a few terms in the Fourier expansion for different 
values of the available energy in the scattering process. 

In this section a new analysis of the gluon Green's 
function as obtained from the NLL BFKL kernel has been 
presented. The method of solution is based on the Monte 
Carlo integration of the phase space of different partonic 
configurations in the multi-Regge and quasi-multi-Regge 
kinematics. This method has many advantages with re- 
spect to previous analysis of the same problem. It allows 
for a reliable study of angular dependences in a straight- 
forward manner, the multiplicities in the evolution are un- 
der control, and it provides an exact solution even with 
running coupling terms which break the scale invariance 
in the kernel. Many other studies are on their way using 
this procedure, as for example, deep inelastic scattering, 
the non-forward case and the matching of this solution to 
different impact factors for the final calculation of cross- 
sections at NLL where the BFKL approach will be relevant 
at present and planned colliders. 



2.5 Resummation at small x 

Main author A. Stasto 

The large magnitude of the NLLx correction in the 
high energy limit, as well as the instabilities associated 
with it, motivate the study of the resummation procedure 
in the limit of small x. In particular it has been observed 
that, by taking into account coUinear limits correctly in 
the NLLx equation, as it is required by the DGLAP dy- 
namics, stabilizes the high energy expansion. To under- 
stand this in more detail let us recall the structure of the 
LLx BFKL equation in the Mellin space where the Mellin 
variable 7 is conjugated to the logarithm of the transverse 
momentum Ink^/A^ 



x(°)(7) = 2^(1) - ^(7) -^-(1- 7) 
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Fig. 9. Projections on conformal spins of the N=4 SUSY NLL 
BFKL Green's function. 



where in the pole expansion of the kernel eigenvalue 
we have retained only leading collinear and anticoUinear 
poles. These correspond exactly to the DGLAP strong or- 
dering of transverse momenta along the gluon ladder. In 
the NLLx case the eigenvalue function takes on a compli- 
cated functional form which in the collinear limit is 



X^^\l) 



AM , ^i(O) 
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7 1 
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■7 
(26) 
with Ai(0) = -11/12. Note the negative sign of the NLLx 
contribution. It turns out that the collinear approximation 
above reproduces the exact result within ~ 7% of accu- 
racy. The terms proportional to ^i(O) are related to the 
non-singular in x part of the LO DGLAP splitting func- 
tion, whereas the cubic poles come from the energy scale 
choice. The highly singular form of the NLLx correction as 
it is seen from eq.(j211) is the source of the large correction 
and potentially unstable behavior. The resummation pro- 
cedure presented in [15] is based on four key ingredients: 

— Taking into account the full splitting function at LO 
in the DGLAP approximation. 

— Incorporating the energy scale change in the form of 
the kinematical constraint. 
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Fig. 10. cr^ot^ in the framework of BFKL 



— Running of the coupling constant as 

— Subtraction of the double and single poles in order to 
avoid double counting. 

In [68] a procedure based on the numerical solution of 
the BFKL equation in momentum space was presented. 
It takes into account all of the above-mentioned ingre- 
dients and yields stable result for the intercept and the 
gluon Green's function. Furthermore, the procedure for 
extraction the resummed splitting function was also pre- 
sented, which is more relevant for application to the deep 
inelastic scattering processes such as measured at HERA. 
In Fig l2.5l we show the resummed splitting function ob- 
tained in the resummed scheme [68], together with the 
renormalization scale variation and the singular in x part 
of the NNLO DGLAP splitting function. The characteris- 
tic feature of the resummed splitting function Pgg is the 
strong preasymptotic behavior at intermediate values of 
X ~ 10"'^ — 10^'' which manifests itself in the dip of the 
splitting function, only later followed by the increase at 
very small x. Also interesting is the fact that the small 
X part of the NNLO DGLAP splitting function matches 
nearly exactly with the initial decrease of the resummed 
splitting function. The existence of the dip rather than 
an increase at values of x ~ 10~^ can have an interesting 
impact on the phenomenology. 



2.6 The NLO 7* impact factor 

Main author A. Kyrieleis 

One of the most attractive observables to test the 
BFKL approach is the total cross section for 7*7* scat- 
tering. To calculate this observable in the framework of 
NLO BFKL the 7* impact factor (^) at NLO is needed in 
addition to the universal BFKL Green function (G), see 
FigHni 



If the NLO BFKL equation is solved in the momentum 
space the numerical value of the 7* impact factor has to 
be known as a function of the Reggeon momentum and of 
the energy scale. 

Besides this, the NLO 7* impact factor allows to ap- 
proach the resummation of the next-to-leading logs(l/x) 



in the quark anomalous dimensions. It also provides the 
full information necessary to investigate the color dipole 
picture at NLO which, at LO, is one of the important in- 
gredients to the QCD evolution based upon the Balitsky- 
Kovchegov equation (see section 5 below). At the first 
small-x workshop [22] first steps in the calculation of this 
impact factor have been presented. 

The virtual and the real corrections of the 7* impact 
factor are calculated from the photon-Rcggeon vertices 
for qq and qqg production, respectively. Both vertices are 
known [88-91]. What remains to complete the calculation 
of the NLO photon impact factor after the infrared di- 
vergences of the virtual and of the real parts have been 
combined [91] are the integrations over the qq and qqg 
phase space, respectively. 

Recently, the phase space integration in the real cor- 
rections have been performed for the case of longitudinal 
photon polarization, [92]. The integration over the trans- 
verse momenta have been carried out analytically. To this 
end the Feynman diagrams were treated separately giving 
rise to additional divergences that have been regularized. 
As the result, a convergent Feynman parameter integral 
has been obtained for each Feynman diagram (or small 
groups of them). These results can serve as a starting 
point for further analytical investigations, in particular 
because the Mellin transform of the real corrections w.r.t 
the Reggeon momentum can be easily obtained. 

The remaining integrations in the real corrections (lon- 
gitudinal 7* polarization) have been carried out numeri- 
cally [92]. The result is a function #'"*=°' of two dimen- 
sionless (scaled by the photon virtuality) variables: the 
Reggeon momentum r^ and the energy scale sq. A physi- 
cal scattering amplitude (e.g. for the 7*7* scattering pro- 
cess) involving the BFKL Green's function and the impact 
factors has to be invariant under changes of sq. The sq de- 
pendence of the 7* impact factor therefore represents an 
important issue, sq enters the NLO 7* impact factor as a 
cutoff to exclude that region of the qqg phase space where 
the gluon is separated in rapidity from the qq pair (LLA) . 
The virtual corrections are therefore independent of sq 
and the integration of the real corrections alone already 
allows to study the sq dependence of the NLO 7* impact 
factor. Let us define, as part of the full NLO impact factor: 
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where (7^^(°) denotes the LO 7* impact factor and g^ = 
Airas- Choosing Q^ — 15 GeV^ for the photon virtuahty 
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Fig. 11. <?' at different different values of sq 



leads as{Q^) = 0.18 01 g — 1.5. Fig ^] compares <!>' to the 
LO impact factor as function of r^ at different values of sq. 
The real corrections are negative and rather large. More 
important, #' becomes, in absolute terms, more significant 
for smaller values of sq. This implies that the 7* impact 
factor tends to become smaller with decreasing sq. Since a 
decrease of sq in the energy dependence ( — )'^ will enhance 
the scattering amplitude, the combined sq dependence of 
the impact factors and the BFKL Green's function has to 
compensate this growth. The result for the sq behavior of 
the 7* impact factor is therefore, at least, consistent with 
the general expectation. To check the sq (in) dependence 
of the full scattering amplitude and to compute cr''' ''' , at 
least for longitudinal 7* polarization, the phase space inte- 
gration in the virtual corrections is the only piece missing. 



3 Applications of /cx-factorization 

In coUinear factorisation the transverse momenta of the 
incoming partons are neglected whereas they are included 
in /c^-factorization if the same order in as of the calcu- 
lation is considered. Thus in coUinear factorsiation these 
transverse momentum effects come in as a next-to-leading 
order level. 

In the following sections we discuss some applications 
of A:j^-factorization to describe heavy quark production in 
pp collisions. 



3.1 Heavy quark production at the Tevatron 

Main author N. Zotov 

Heavy quark production in hard collisions of hadrons 
has been considered as a clear test of perturbative QCD. 



Such processes provide also some of the most important 
backgrounds to new physics phenomena at high energies. 
Bottom production at the Tevatron in the k±- 
factorization approach was considered earlier in [10,28,93- 
98]. Here we use the /c^-factorization approach for a more 
detailed analysis of the experimental data [36,38,99-101]. 
The analysis also covers the azimuthal correlations be- 
tween b and b quarks and their decay muons. Some of these 
results have been presented earlier in Refs. [98, 102-107] 
(see also [22,23]). 



3.2 Theoretical framework 

In the fcj^-factorization approach, the differential cross sec- 
tion for inclusive heavy quark production may be written 
as (see [108]) 



da{pp^QQX)^ 
1 



[9*9* - Q 

X dyi dy2 dPar aQiT d^.2T 1^ 1^ 



167r(a;i X2 s)' 
X E \M\W{9*9* - QQ) 



d4i_ d(h d(f>Q 
1-K 1-K 2tt 



(27) 



where A{xi,o^rp,^^) and ^(x2,q2T,M^) are unintegrated 
gluon distributions in the proton, qiT, q2T, P2T and 0i, 
4>2, 4>Q are transverse momenta and azimuthal angles of 
the initial BFKL gluons and final heavy quark respec- 
tively, 2/1 and j/2 are the rapidities of heavy quarks in 
the pp center of mass frame. ^ l-^lsHA(ff*5* ~^ QQ) i^ 
the off mass shell matrix element, where the symbol ^ 
in 127(1 indicates an averaging over initial and a summa- 
tion over the final polarization states. The expression for 
X] l^'-^liHA(5*5* ^ QQ) coincides with the one presented 
in [11]. 

In the numerical analysis, we have used the KMS pa- 
rameterization [52] for the /cx-dependent gluon density. It 
was obtained from a unified BFKL and DGLAP descrip- 
tion of F2 data and includes the so called consistency con- 
straint [17]. The consistency constraint introduces a large 
correction to the LO BFKL equation; about 70% of the 
full NLO corrections to the BFKL exponent A are effec- 
tively included in this constraint, as is shown in [17, 109]. 



3.3 Numerical results 

In this section we present the numerical results of our 
calculations and compare them with i3-meson production 
at DO [36,101], CDF [38,99,100] and UAl [110]. 

Besides the choice of the unintegrated gluon distribu- 
tion, the results depend on the bottom quark mass, the 
factorization scale p? and the b quark fragmentation func- 
tion. As an example, Ref. [Ill] used a special choice of the 
6-quark fragmentation function, as a way to increase the 
B meson cross section in the observable range of trans- 
verse momenta. In the present paper we convert b quarks 
into B mesons using the standard Peterson fragmentation 
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Fig. 12. The b quark transverse momentum distribution (in- 
tegrated from Pt i-nin) &t Tevatron conditions presented in the 
form of integrated cross sections. The curves correspond to 
the fcx-factorization results with the KMS unintegrated gluon 
distribution. Experimental data are from UAl [110] (Fig. a)), 
DO [36] (Fig. b)), and CDF [38,99] (Fig. c)). 



function [112] with e = 0.006. Regarding the other param- 
eters, we use rub = 4.75 GeV and /i^ = qf. as in [10, 95]. 

The results of the calculations are shown in Figs. I12t 
lltil Fig. 1121 displays the b quark transverse momentum 
distribution at Tevatron conditions presented in the form 
of integrated cross sections. The following cuts were ap- 
plied: (a) Jyi] < 1.5, Jysl < 1-5, ^ = 630 GeV; (b) 
\yi\ < 1, Vi = ISOOGeV; and (c) \yi\ < 1, Jyzl < 1, 
^/s — ISOOGeV. One can see reasonable agreement with 
the experimental data. 

Fig. 1131 shows the prediction for the B meson px spec- 
trum at Vs = 1800 GeV compared to the CDF data [38] 
within the experimental cuts \y\ < 1, where also a fair 




Fig. 13. Theoretical predictions for the B meson pr spec- 
trum compared to the CDF [100] data. Curve is the same as 
in Fig. ini 



agreement is found between results obtained in the k±- 
factorization approach and experimental data. 

The DO data include also muons originating from 
the semileptonic decays of _B-mesons. To produce muons 
from B mesons in theoretical calculations, we simulate 
their semileptonic decay according to the standard elec- 
troweak theory. In Fig. 1141 we show the rapidity distribu- 
tion da/dly'^l for decay muons with pt^ > 5 GeV. 

Fig. ll5l shows the leading muon px spectrum for bb pro- 
duction events compared to the DO data. The cuts applied 
to both muons are given by 4 < p^ < 25 GeV, jTy^J < 0.8 
and 6 < m'^^ < 35 GeV. The leading muon in the event 
is defined as the muon with largest p^-value. In all the 
above cases a rather good description of the experimental 
measurements is achieved. 

It has been pointed out that investigations of bb corre- 
lations, such as the azimuthal opening angle between b and 
b quarks (or between their decay muons) , allow additional 
details of the b quark production to be tested, since these 
quantities are sensitive to the relative contributions of the 
different production mechanisms [10,28,93-95,97]. In the 
collinear approach at LO the gluon-gluon fusion mecha- 
nism gives simply a delta function, S{A(j)'''' — tt), for the 
distribution in the azimuthal angle difference Acf)''^. In the 
fcj^-factorization approach the non-vanishing initial gluon 
transverse momenta, qiT and q2T, implies that this back- 
to-back quark production kinematics is modified. In the 
collinear approximation this effect can only be achieved if 
NLO contributions are included. 

The differential bb cross section da/dAcf)'^^ is shown 
in Fig. El (from [108]). The following cuts were applied 
to both muons: 4 < p^ < 25 GeV, \t]''\ < 0.8 and 6 < 
m'^^ < 35 GeV. We note a significant deviation from the 
pure back-to-back production, corresponding to Aip'^^ « 
TT. There is good agreement between the KMS prediction 
and the experimental data, which shows that for these 
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Fig. 14. The cross section for muons from B meson decay as a 
function of rapidity compared to the DO data [101]. The curves 
are the same as in Fig. 1121 
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Fig. 15. Predictions for the leading muon pr spectrum in the 
bb production events compared to the DO data [36]. The curve 
is the same as in Fig. 1121 



correlations the fc_L-factorization scheme with LO matrix 
elements very well reproduces the NLO effects due to the 
gluon evolution. 




Fig. 16. Azimuthal muon-muon correlations at Tevatron con- 
ditions. The curve is the same as in Fig. 1121 Experimental data 
are from the DO collaboration [36]. 



3.4 Quarkonium production 

Main author S. Baranov 

The fc_L-factorization approach has rather successfully 
described the production of open charm and beauty, as 
discussed in the previous section, but also hadroproduc- 
tion of heavy quarkonium states, J/V', Xc a-nd T mesons, 
at the Tevatron are well described [95, 113-115]. In many 
cases, however, the data can also be described within 
the usual coUinear parton model, if the relevant next-to- 
leading order QCD corrections are taken into account, or 
if the so called color-octet mechanism is included. 

In this context, the theoretical predictions on J /'ip spin 
alignment made in Ref. [116] are of particular interest, as 
the coUinear and fc_L -factorization approaches show quali- 
tatively different behavior. Note that the /s^-factorization 
approach provides the only known (up to date) explana- 
tion of the J /ip polarization phenomena observed at the 
Tevatron [117] and at HERA [105]. 

It would be interesting and important to find other 
examples, where the difference between the coUinear and 
noncoUinear approaches would be manifested in a clear 
and unambiguous way. In this section we suggest such 
a process. We analyze the production of P-wave quarko- 
nium states (namely the Xc and Xb mesons) in high energy 
hadronic collisions and demonstrate the dramatic differ- 
ence between the different theoretical calculations. 

Naively one could expect a difference from the fact 
that the production of xi states in the 2^-1 gluon-gluon 
fusion process is forbidden, if the initial gluons are on 
shell, but is allowed if the gluons are off shell. However, 
the real situation is complicated by the necessity to take 
into account also the 2^-2 processes. The results of our 
analysis are presented in the next subsection. 

We begin our discussion with showing the predictions 
of the coUinear parton model for the production of P- 
wave charmonia at Tevatron conditions. The color-singlet 
production scheme refers to the 2 — > 2 gluon-gluon fusion 
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subprocess 



9 + g ^X + 9- 



(28) 



(It would be inadequate to rely upon the 2^-1 sub- 
process g + g — > X in this case, because the final state 
particle would then be produced with zero transverse mo- 
mentum, and thus could not be detected experimentally.) 
The computational technique is explained in detail else- 
where [118-120]. 

For the sake of definiteness, we only present the pa- 
rameter setting used in our calculations. Throughout the 
paper we use the LO GRV set [121] for gluon densities 
in the proton, and the value for the Xc wave function, 
|7^'^J0)|2 = 0.075 GeV^ taken from the potential model 
of Ref. [122]. The renormalization scale in the strong 
coupling constant as(/z|j./yl^) is set to /x^ — rn'^ + Pr x 
with yl=200 MeV. The integration over the final state 
phase space is restricted to the pseudorapidity interval 
—0.6 < ri{xc) < 0.6, in accord with the experimental cuts 
used by the CDF collaboration [123-128]. 

Since in the coUincar formalism the predictions based 
on the color-singlet mechanism alone are known to be in- 
consistent with the data [123-128], the theory has to be 
amplified with the so called color-octet contribution, as 
it is commonly assumed in the literature [120]. Unlike the 
predictions of the color-singlet model, the size of the color- 
octet matrix elements are not calculable within the theory. 
Therefore, the corresponding numerical results are always 
shown with arbitrary normalizing factors (just chosen to 
fit the experimental data when possible). 

The numerical predictions of the coUinear parton 
model are summarized in Fig. 1171 fupper panel). At rela- 
tively low transverse momenta, the production of Xc states 
is dominated by the color singlet mechanism. The differ- 
ential cross section da/dpr diverges when pT — > for X2 
states (dashed histogram), while it remains finite for xi 
states (solid histogram). The production of xi states at 
zero pt is suppressed (in accord with the Landau- Yang 
theorem) , because in the limit of very soft final state glu- 
ons the 2^2 gluon-gluon process degenerates into the 
2^1 process. The shape of the xo spectrum is similar to 
that of X2 (up to an overall normalizing factor), and this 
spectrum is not shown in the figure. 

The production of Xc mesons at high px is dominated 
by the color-octet contribution, which mainly comes from 
the 'gluon fragmentation' diagrams. Here, the perturba- 
tive production of ^Si color octet states. 



9 + 9 



3c8 



'St 



(29) 



is followed by a nonperturbative emission of soft gluons, 
which results in the formation of physical color singlet Xc 
mesons: 

3c8 , 3pl 



^5f 



J 



ng. 



(30) 



As the co-produced gluons in cq. (|30|l are assumed to 
be soft, the momentum distribution of Xc mesons is taken 
identical to that of the color-octet '^Si state in eq. (|29|l . 
The nonperturbative matrix elements responsible for the 
process eq. (|30|l are related to the fictitious color-octet 
wave functions, which are used in calculations based on eq. 



(|29|l in place of the ordinary color-singlet wave function: 
<0lO8l0>=(9/27^)l7^8(0)l2. 

It should be noted that the fragmentation of an almost 
on-shell transversely polarized gluon into a xi state via the 
emission of a single additional gluon, g -^ ^Sf ^ Xi -I- <7, 
is suppressed in accord with Landau- Yang theorem. In 
terms of the nonrelativistic approximation, it is equiva- 
lent to say that the formally leading color-electric dipole 
transitions are forbidden, and one must go to nonlead- 
ing higher multipoles. As the degree of this suppression is 
not calculable within the color-octet model on its own, we 
rather arbitrarily set the suppression factor to 1/20, which 
corresponds to potential model expectations for the aver- 
age value of u^ . 

We now proceed with showing the results obtained in 
the /c^-factorization approach. In this case the production 
of charmonium Xc states can be successfully described 
within the color-singlet model alone [117], or with only 
a minor admixture of color-octet contributions [95] . The 
consideration is based on the 2^1 partonic subprocess 



9 



X, 



(31) 



which represents the true leading order in perturbation 
theory. The nonzero transverse momentum of the final 
state meson comes from the momenta of the initial gluons. 
The computational technique, which we are using here, is 
identical to the one described in detail in Ref. [117]^. 

In order to estimate the degree of theoretical uncer- 
tainty connected with the choice of unintegrated gluon 
density, we also use the prescription proposed in [9]. In 
this approach, the unintegrated gluon density is derived 
from the ordinary density G(x, q^) by differentiating it 
with respect to q^ and setting q^ = k\. Among the differ- 
ent parameterizations available on the present-day theo- 
retical market, this approach shows the largest difference 
with Bliimlein's density [129]. Thus, these two gluon den- 
sities can represent a theoretical uncertainty band. 

The numerical results are exhibited in Fig. 1171 fmiddle 
panel). In contrast with the coUinear parton model, the 
differential cross sections are no longer divergent, even 
at very low pT values. This property emerges from the 
fact that the relevant 2^1 matrix elements are al- 
ways finite. One can see that the production of the xi 
state (solid histogram) at low px is strongly suppressed 
(in comparison with the xo ^-nd X2 states, short and long 
dashed histograms) because the initial gluons are almost 
on-shell. The suppression goes away at higher pT, as the 
off-shellness of the initial gluons becomes larger. 

In Fig. 1171 flower panel) we compare the predictions of 
the coUinear and fcj^-factorization approaches by showing 
the ratio of the differential cross sections da{xci)/dpT and 
da{xc2)/dpT plotted as a function of pt- As long as the 
ratio of the nonperturbative color-octet matrix elements, 
0{^Sf -^ xi)/0{^Sf ^ X2), is unknown, the predictions 
of the coUinear parton model are very uncertain. The dif- 
ferent dotted curves in Fig. El from top to bottom corre- 
spond to the color-octet X1/X2 suppression factor set to 1, 

^ We use the FORTRAN code developed in [117]. This code 
is public and is available from the author on request. 
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Fig. 17. Theoretical predictions for the production of Xc 
mesons at Tevatron conditions. Upper: Predictions of the 
coUinear parton model. Solid histogram, xi production via 
color-singlet mechanism; dashed histogram, X2 production via 
color-singlet mechanism; the lower and the upper dotted his- 
tograms, Xi s-nd X2 production via color-octet mechanism, 
respectively. Middle: Predictions of the fcx-factorization ap- 
proach. Solid histograms, xi production; thin and thick dashed 
histograms, xo and X2 production, respectively. The upper and 
the lower histograms of each type correspond to the gluon den- 
sities of Refs. [129] and [9]. Only the color singlet mechanism 
is assumed in all cases. Lower: Predictions on the ratio of 
the differential cross sections dcr(xi)/rfcr(x2)- Solid histograms, 
fcx-factorization approach with gluon densities of Refs. [129] 
and [9]; dashed histogram, coUinear parton model, color singlet 
contribution only; dotted histograms, coUinear parton model 
with both singlet and octet production mechanisms taken into 
account. The different curves from top to bottom correspond 
to the color-octet Xi/X2 suppression factor set to 1, 0.3, 0.1 
and 0.03, respectively. 



Fig. 18. Theoretical predictions on the production of Xb- 
notations are the same as in Fig. 1171 



The 



0.3, 0.1, and 0.03, respectively. The band between the two 
lowest histograms may be considered as the most realistic 
case. The predictions of the coUinear and fcj^-factorization 
approaches clearly differ from each other in their absolute 
values, and show just the opposite trend in the experimen- 
tally accessible region {px > 5 GeV). 

We conclude our discussion with showing the predic- 
tions for the bottomonium states. The calculations are 
performed with the parameter setting given above, and 
with the value of the Xb wave function set equal to 
|7^^^(0)|2 = 1.4 GeV^ [130]. The integration over the final 
state phase space is now restricted to the pseudorapid- 
ity interval —0.4 < ri{xb) < 0.4, in accord with the CDF 
experimental cuts [123-128]. 

Our numerical results are displayed in Fig. 1181 The 
qualitative features of the differential cross sections are 
similar to the ones, which we have seen in the case of 
charmonium. It is worth recalling that the production of 
T mesons has been already measured by the CDF col- 
laboration [123-128] at Pt values close to zero. Although 
the Pt dependence of the direct {pp — > TX) and indirect 
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{pp -^ XbX ^ "fl^) contributions have not been stud- 
ied separately, the net result seems to be at odds with 
collinear calculations. In fact, the predicted magnitude of 
the indirect contribution coming from the decays of X62 
states at pt < 2 GeV exceeds the total measured T pro- 
duction rate in this region. In contrast the measured differ- 
ential cross section da{T)/dpT decreases with decreasing 
PTj in perfect agreement with the fc_L-factorization predic- 
tions [117]. 

In summary, one major difference between the collinear 
and the fci^-factorization approaches is connected with the 
behavior of the differential cross section da{x2) / dpr at 
low transverse momenta. This quantity remains finite in 
the fc^-factorization approach, while it diverges in the 
collinear parton model when pT goes to zero. The lat- 
ter prediction seems to be not supported by the available 
experimental data on the bottomonium production at the 
Tevatron. 

Another well pronounced difference refers to the ratio 
between the production rates da{xi)/d<j{x2)- The under- 
lying physics is connected with the off-shellness of the glu- 
ons. In the collinear parton model the relative suppression 
of xi states becomes stronger with increasing pT because 
of the increasing role of the color-octet contribution. In 
this approach the leading-order fragmentation of an on- 
shell transversely polarized gluon into a vector meson is 
forbidden. In contrast with that, in the A:j_-factorization 
approach the increase in the final state px is only due 
to the increasing transverse momenta (and correspond- 
ing virtualities) of the initial gluons, and consequently the 
suppression motivated by the Landau- Yang theorem be- 
comes weaker at large pr ■ 

In conclusion we see that quarkonium production can 
be regarded as a direct probe of the gluon virtuality, and 
provides a direct test of the need for a noncoUinear parton 
evolution. Our results seem especially promising in view 
of the fact that the difference between the two theoretical 
approaches is clearly pronounced at conditions accessible 
for direct experimental measurements. 



4 BFKL dynamics in jet-physics 

Main author G. Marchesini 

It has been generally taught that QCD dynamics in 
high-energy scattering and in jet-physics are quite differ- 
ent. However it has been recently shown [131] that classes 
of jet observables satisfy equations formally similar to the 
ones for the high-energy ^'-matrix. The jet-physics observ- 
able here discussed are the heavy quark-antiquark multi- 
plicity (in certain phase-space region) and the distribution 
in the energy emitted away from jets. They satisfy equa- 
tions formally similar to BFKL and Kovchegov equations 
respectively. One may expect that by exploiting such a 
formal similarity will bring new insights in both fields. 

The common key feature shared by the observables 
in these two cases is that enhanced logarithms come only 
from infrared singularities (no collinear singularities). The 
differences between the two cases is in the relevant phase 



space for multi soft-gluon ensemble. For the S'-matrix all 
transverse momenta of intermediate soft gluons are of 
comparable order (no collinear singularities in transverse 
momenta). For the considered jet-observables all angles of 
emitted soft gluons are of comparable order (no collinear 
singularities in emission angles). 

We discuss first the QQ (heavy quark-antiquark) mul- 
tiplicity in the phase-space region where collinear singular- 
ities cancel and then the distribution in the energy emitted 
away from jets. 



4.1 QQ-multiplicity and BFKL equation 

The standard multiplicity in hard events has both 
collinear and infrared enhanced logarithms which are re- 
summed by the well known expression [132, 133]. 



dkt 






Ncas 



(32) 



The QQ-multiplicity introduced and studied in [131] is, 
due to the peculiar phase space region chosen, without 
collinear singularities. In e^e" with center of mass en- 
ergy Q one considers the emission of a QQ system of mass 
M. and momentum k. In the calculation one takes: small 
velocity v = \k\/Ek so that there are no collinear singu- 
larities; Q ^ A^ so that perturbativc coefficients are en- 
hanced by powers oilnQ/M; and studies the process near 
threshold. In this region, the leading logarithmic contribu- 
tions (a" In" Q/A4) are obtained by considering soft sec- 
ondary gluons qi,- ■ ■ Qn emitted off pp, the primary quark- 
antiquark. The QQ system originates from the decay of 
one of these soft gluons, actually the softest one, we de- 
note by fc. 



pp + qi.. .qnk, 



(33) 



As shown in [131], to leading logarithmic order, the QQ- 
multiplicity distribution factorizes into the inclusive dis- 
tribution / for the emission of the soft off-shell gluon of 
mass M and momentum \k\ and the distribution for its 
successive decay into the QQ system 



EkdN 



dM^d\k\ Stt^M^^ 



aiCp IM^ -AW^M^ + 2A'P 



M2 



M"^ 



/, (34) 



where M is the heavy quark mass. The Born distribution 
is 



47r 



Wab{k) 



(ab) 



{ak){kh) 



{I -cos Bab) 



(1 - vcos9ak){'^ - VCOsOkb)'' 



(35) 



with Wab{k) the (angular part of the) distribution for the 
off soft gluon emitted off the a6-dipole (for e+e^ in center 
of mass 9ab — tt). For v <1 the Born contribution is finite. 
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For Q ^ A4, secondary radiation contributes. Since 
the Born contribution is regular, only soft logarithms 
(a" In" Q/A4) are generated which need to be resummed 
by recurrence relation. To understand the structure of the 
resulting equation and appreciate the similarity with the 
BFKL equation we consider the first non trivial contribu- 
tion in which, besides the off-shell soft gluon k, there is an 
additional massless soft gluon either emitted or virtual. 

The real emission contribution is given by 



Wabik;q) 



{ab) 



(ab) 



{aq){qk){kb) {ak){kq){qb) 

=&{q-k) Wab{q) ■ [Waq{k) + Wqb{k)] 
+0{k-q) Wab{k) ■ [wak{q) + Wkb{q)] 



(36) 



where, for massless q, 
(ab) 



Wab{q) 



1 — cos t 



{aq){qb) {I ~ cos 6 aq) {I -cob 6 qb) 



(37) 



The corresponding virtual correction is obtained by inte- 
grating over the massless momentum q in the expression 
(softest gluon emitted off external legs) 



^Ibif^'^l) = - 0{q-k) Wab{q) ■ Wab{k) 

- 0{k-q) Wab{k) ■ [wakiq) + Wkb{q)] 

By summing the two contributions one finds 
w^+''{k;q)=0{q-k)wab{q) 

X [Waqik)+Wqb{k)-Wab{k)] , 



(38) 



(39) 



which shows that k is the softest gluon. From this we 
derive the first iterative structure giving I^^' in terms of 
the Born contribution (I35II 



l'-^Hpab,T) = / —^as{qt] 
Jm It 



dQq 

47r 



Wab{q) 



I^'Hpaq)+I^"\pqb)-I^'Hpab) 
1 — COS0,-, 



(40) 



P^3 - 



with 



f'^ dqt _ , 
JM It 



2N? 



UNc - 2nf 



■In 



lnQ/A \ 
InM/AJ 



(41) 



Here the running coupling in qt is restored so r is given 
by an expansion in as{Q)\nQ/Ai. The measure in H4()|l 
is the branching distribution for a massless soft gluon q 
emitted off the a6-dipole. One generalizes this branching 
structure as successive dipole emission of softer and softer 
gluons and one deduces [131] 



drl{pab, t) 



I 



-^Wab{q) I{paq,T)+I{pqb,T)~I{pab,T) 



(42) 



This recurrence structure is very similar to the one ob- 
tained in the dipole formulation of the BFKL equation 



[134-136]. The fundamental difference is that here the in- 
clusive distribution / depends on the angular variable p 
(with the limitation p < 1), while in the high energy scat- 
tering one deals with the S-matrix as a function of the 
impact parameter b (which is not bounded). 

The similarity with the BFKL equation can be made 
even more evident if one performs the azimuthal integra- 
tion. One obtains [131] 



drl{p, r)^f^ {v-'livP, r) - I{p, r)) 
Jo ^-V 

+ f7^{nV-'p,r)~Iip,r)). 



(43) 



The lower limit 77 > p in the second integral ensures that 
the argument of I{p/ri, t) remains within the physical re- 
gion p/77 < 1 . The presence of this lower bound is the only 
formal difference with respect to the BFKL equation for 
the high energy elastic amplitude T in the impact param- 
eter representation 

drT{p, r)^f^ {v-'nrjp, r) - T{p, r)) 

-'\ "^ (44) 

+ / -^{nr'p,T)-T{p,T)). 
Jo ^-V 

Here p = 6^ is the square of the impact parameter and 
T = oisY with Y the rapidity with the QCD coupling 
fixed. We discuss now the differences in the two solutions. 
Recall first the solution for the high-energy scattering 
case. Since b has no infrared bound we change the variable 



b^ = e- 



— cx) < a; < 00 . (45) 

The BFKL equation 1)44(1 satisfies translation invariance 
and the area conservation law 

r°° 1 

dr dxe2^r(e-^r)e-''i"2^ = 0. 

J —00 

This allows us to obtain the solution and then its asymp- 
totic behavior (using D = 28^(3) = 33.6576 . . .) 



T{b,T) 



,27r 



T(fc)e(*'=-2)^e 



.f(0) 



g41n2r g-jXg-^Ej^ 



(46) 



^/2ttDi 



with T determined by the initial condition and x(fc) = 
2ip{l)—ip{^+k)—4>{^—k) the BFKL characteristic function. 
In the QQ-multiplicity case, the crucial difference is 
that the angular variable p is bounded. Introducing the x- 
variable as in (|45|l one observes that translation invariance 
is lost and, instead of area conservation, one has absorp- 
tion 

1 — cos^ _ 
p= =e , < a; < 00 , 

/■°° 1 

Jo 
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Fig. 19. Plot of 0(x,r) = e"2^ e-'""2^J(e-^, r) solution of 
14311 with initial condition I{p, 0) — ^p. 



The exact solution of (|^ was obtained in [137] 



IiP,r) 



dki{k)P^i 

2 

(a;+Xo)e'*''^2r 



■+ik 



P] pxWr 



(47) 



-V2ttDi 



with Pa{z) the Legendre function (well known in Regge 
theory) and / given by initial condition. 

From H47|) and from the upper plot of Fig. ^1 one 
has that the inclusive distribution vanishes at the non- 
physical point X = —xq which is slowly varying with r. 
The asymptotic shape is developed already at relatively 
small T. At X = 0, corresponding to the physical value 
p = 1 for e^e~ in center of mass, the function (J){x,t) 
is decreasing, however, thanks to the e'*'"^'^ factor the 
inclusive distribution I{p = 1, r) is increasing as shown in 
the lower plot of Fig. ^1 



4.2 Away-from-jet energy flow in e+e 

Consider in e+e^ annihilation the distribution in the en- 
ergy emitted outside a cone around the jets, -Eout: 



out 




^e+e- (^out) = 



=E 



dan 



0\E, 



-E%» 



contains only single logarithms 
coming from soft singularities so 



This is the simplest (in principle) observable involving 
non-global single logarithms which were (re) discovered by 
Mrinal Dasgupta and Gavin Salam [138-141]. These enter 
all jet-shape observables which involve only a part of phase 
space and therefore are present in a number of distribu- 
tions such as: Sterman- Weinberg distribution (energy in a 
cone); photon isolation; away from jet radiation; rapidity 
cuts in hadron-hadron (e.g. pedestal); DIS jet in current 
hemisphere. As for the observable previously discussed, 
these non-global logs originate from multiple soft gluon 
emissions at large angles (i.e. not in coUinear configura- 
tion). 

^e+e- (£'out) 

«ln"0/£;out) 
that ct„/(Tt can be taken as the distribution in the 
number of soft gluons emitted off the primary pp quark- 
antiquark pair which is known [142] in the large Nc limit. 
Z'e+e-(£^out) was first studied [138-141] numerically by a 
Monte Carlo method and then studied [143] analytically 
by deriving the following evolution equation 

dr^ab = —{drRab) Sab+ / — ; Wab{q) [^aq ' ^qb ~ -S'ab] , 

Jin 47r 

[ dn„ 

Rab = T / — Wab[q) , 

(48) 

where r is the single logarithmic variable previously intro- 
duced (I41|) . As before, to set up a recurrence relation, one 
needs to generalize the problem by introducing distribu- 
tion Sab — ^ab{Eout) for the emission off a6-dipole form- 
ing an angle 9ab- The physical distribution Z'g+g- (i?out) 
for e+e~ in the center of mass is obtained by setting 

Oab = TT. 

As shown in (|48|l . the dipole directions a and b are 
always inside the jet region (g in the integral is bounded 
inside the jet region). If a,b are in opposite semicones, 
then either a, q oi q,b are in the same semicone. There are 
many properties of this jet-physics equation (see [138-141] 
and [143]). What concerns us here as far as the connection 
with high-energy physics is the case in which a, b are in the 
same semicone and we consider a very close to b. In the 
small angle limit we introduce the 2-dimensional variable 
6 for the a6-dipole {Sab ^^ ^(^))- For small 6 we can 
neglect the linear term {Rab ^ S^) so that the evolution 
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equation H48|l becomes 



i(^-(^y (49) 

s{T,e')E{T,e-e')-E{T,e)\, 



with 0' ranging in the full plane. The initial condition is 
17(0,0) = 1. This equation is formally the same as the 
Kovchegov equation [144] for the S-matrix 



drS{T,b) 



'2h' 



d^b 



2- b'\b-b'r (50) 

'S{T,b')SiT,b-b')-S{T,b)], 



where b is the impact parameter ranging in the full plane 
and T = oigY as before. Here the initial condition is 1 — 
5(0, 6) ^ ttg corresponding to the two gluon exchange. 

The asymptotic properties of the solutions are well 
known. Both solutions undergo well known saturation for 
the variable 9^ or b^ larger than a critical value with 
asymptotic behavior e~'^'^ with c ~ 4.88 • • • determined 
from the BFKL characteristic function. Beyond such a 
critical value the solution decreases in r as a Gaussian, 

The difference in the initial condition makes a dif- 
ference in the way the saturation regime is asymptoti- 
cally reached in the two cases. In the high-energy case 
(1 — 5'(0,b) ~ al) the saturation regime is reached after 
a critical time Tc ^ lnQ;7^/41n2. In the jet-physics case 
(17(0, 6) = 1) there is not a critical t and the solution 
goes without impediment into the saturation regime. 

In addition to the different initial conditions, an im- 
portant difference is that the variables in H48|) are angular 
variables ranging in compact regions. We have seen in the 
previous analysis that even at small angle it is not fully 
correct to neglect the compactness affecting the integra- 
tion limits. This question will be further studied [145]. 

4.3 Physics differences 

The basis for the two classes of equations, (|43|I . H48() in 
jet-physics and (|44|l . 150|l in high-energy scattering, is of 
course the (same) multi-soft gluon-distribution. However 
the dominant contributions for the two classes of observ- 
ables {lab, ^ab and T, S) are obtained from very different 
kinematical configurations as we discuss now. 

Jet-physics case: Here all angles 9i of emitted gluons 
are of same order. This is due to the fact that this observ- 



able does not contain coUinear singularities for 6*, 



0. 



Moreover, in the (leading) infrared limit soft gluon ener- 
gies can be taken ordered so that also the emitted trans- 
verse momenta qu are ordered. The ordered variables qu 
enter the argument of the running coupling. The distribu- 
tion lab or Sab are functions of the angular variable 9ab 
(which ranges in a compact region) and r, the logarithmic 
integral of the running coupling in ()41|l . We arc then in- 
terested in the solution for finite 9ab (e.g. 9ab = tt in e 
center of mass) and for r never too large. 



+^- 



High-energy scattering case: Here all intermedi- 
ate soft gluon transverse momenta qu are of same order 
(no singularities for vanishing transverse momentum dif- 
ferences). On the other hand, energy ordering implies in 
this case that intermediate gluon angles 9i are ordered. 
Contrary to the previous case, the running coupling is a 
function of the variables qu which all are of same order. 
Therefore, in first approximation, one can take a^ fixed. 
The high-energy S'-matrix is a function of the impact pa- 
rameter (which has no bound at large b) and t = agY. In 
this case we are then interested in the solution for small 
p (the short distance region) and for t large. 

As discussed in section ITTl the fact that the variable 
p entering the jet-observable ranges in a compact region 
affects the prefactor of the asymptotic behavior and the 
shape of the distribution at finite angles. In the non linear 
case discussed in section ^21 even neglecting compactness 
at small angle, the difference in the initial conditions af- 
fects the ranges in t at which the asymptotic behavior 
(saturation) is developing. 

Concluding, by exploiting similarities and differences 
in the dynamics of high energy scattering and jet-physics 
(with non-global logs) one hopes that new insights in both 
fields could be developed. 



5 Saturation 

Main authors M. Lublinsky and K. Kutak 

A parton evolution equation which attempts to de- 
scribe saturation phenomena was originally proposed by 
Gribov, Levin and Ryskin [9] (GLR equation) in momen- 
tum space and proven in the double log approximation of 
perturbative QCD by Mueller and Qiu [146]. In the lead- 
ing Inl/x approximation it was derived by Balitsky in 
the Wilson Loop Operator Expansion [147]. In the form 
presented later it was obtained by Kovchegov [144] (now 
called the Balitsky-Kovchegov, or BK equation) in the 
color dipole approach [134] to high energy scattering in 
QCD. This equation was also obtained by summation of 
the BFKL pomeron fan diagrams by Braun [148] and most 
recently Bartels, Lipatov, and Vacca [149]. In the frame- 
work of Color Glass Condensate it was obtained by lancu, 
Leonidov and McLerran [150]. 



5.1 Basic facts about the BK equation 

Because the transverse coordinates are unchanged in a 
high energy collision, unitarity constraints are generally 
more easy to take into account in a formalism based 
on the transverse coordinate space representation, and 
several suggestions for how to include saturation effects 
in such a formalism have been proposed. Golec-Biernat 
and Wiisthoff [151] formulated a dipole model, in which 
a virtual photon is treated as a qq or qqg system im- 
pinging on a proton, and this approach has been further 
developed by several authors (see e.g. [152] and [153]). 
Mueller [134, 135, 154] has formulated a dipole cascade 
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model in transverse coordinate space, which reproduces 
the BFKL equation, and in which it is also possible to ac- 
count for multiple sub-collisions. Within this formalism 
Balitsky and Kovchegov [144, 147] have derived a non- 
linear evolution equation (BK equation), which also takes 
into account these saturation effects from multi-pomeron 
exchange and which is the best presently available tool to 
study saturation phenomena at high energies. Contrary to 
many models the BK equation has solid grounds in per- 
turbative QCD. The equation reads 



dJV(xoi,y; 
dy 



2tt 



d'x2- 



^02 ^12 



2 7V(xo2,?;;b- -xia) - N{xoi,y\h) 
-iV(xo2, y; b - -xi2)iV(xi2, y; b - -X02; 



(51) 



The function N{r±, x; b) is the imaginary part of the am- 
plitude for a dipole of size r± elastically scattered at an 
impact parameter b. 

In the equation H51|l. the rapidity y = — Inx. The ul- 
traviolet cutoff p is needed to regularize the integral, but 
it does not appear in physical quantities. We also use the 
large Nc limit (number of colors) value of Cp = N'c/2. 

Eq. I)51|l has a very simple meaning: The dipole of size 

Xqi decays in two dipoles of sizes X12 and X02 with the de- 

2 
cay probability given by the wave function Itf'P — ■^"^■i . 

^02 ^12 
These two dipoles then interact with the target. The non- 
linear term takes into account a simultaneous interaction 
of two produced dipoles with the target. The linear part of 
eq. (|51|l is the LO BFKL equation [7,8], which describes 
the evolution of the multiplicity of the fixed size color 
dipoles with respect to the energy y. For the discussion 
below we introduce a short notation for eq. (|51() : 



dN 
dy 



as Ker (g) {N - N N) . 



(52) 



The BK equation has been studied both analytically 
[155-161] and numerically [148,162-168]. The theoretical 
success associated with the BK equation is based on the 
following facts: 

— The BK equation is based on the correct high en- 
ergy dynamics which is taken into account via the LO 
BFKL evolution kernel. 

— The BK equation restores the s-channel unitarity of 
partial waves (fixed impact parameter) which is badly 
violated by the linear BFKL evolution. 

— The BK equation describes gluon saturation, a phe- 
nomenon expected at high energies. 

— The BK equation resolves the infrared diffusion prob- 
lem associated with the linear BFKL evolution. This 
means that the equation is much more stable with 
respect to possible corrections coming from the non- 
perturbative domain. 

— The BK equation has met with phenomenological suc- 
cesses when confronted against DIS data from HERA 
[162,165,169-175]. 



The BK equation is not exact and has been derived in 
several approximations. 

— The LO BFKL kernel is obtained in the leading soft 
gluon emission approximation and at fixed a^- 

— The large Nc limit is used in order to express the non- 
linear term as a product of two functions N. This limit 
is in the foundation of the color dipole picture. To a 
large extent the large N^ limit is equivalent to a mean 
field theory without dipole correlations. 

— The BK equation assumes no target correlations. Con- 
trary to the large Nc limit, which is a controllable ap- 
proximation within perturbative QCD, the absence of 
target correlations is of pure non-perturbative nature. 
This assunrption is motivated for asymptotically heavy 
nuclei, but it is likely not to be valid for proton or re- 
alistic nucleus targets. 

There are several quite serious theoretical problems 
which need to be resolved in the future. 

— The BK equation is not symmetric with respect to 
target and projectile. While the latter is assumed to 
be small and perturbative, the former is treated as a 
large non-perturbative object. The fan structure of the 
diagrams summed by the BK equation violates the t- 
channel unitarity. The i-channel unitarity is a com- 
pleteness relation in the i-crossing channel. It basi- 
cally reflects a projectile-target symmetry of the Feyn- 
man diagrams. The down-type fan graphs summed by 
the BK equation, obviously violate the symmetry. A 
first step towards restoration of the t-channel unitar- 
ity would be an inclusion of Pomeron loops. 

— Though the BK equation respects the s-channcl uni- 
tarity ^ the exchange of massless gluons implies that 
it violates the Froissart bound for the energy depen- 
dence of the total cross section. In order to respect the 
Froissart bound, gluon saturation and confinement are 
needed. On one hand, the BK equation provides gluon 
saturation at fixed and large impact parameters. On 
the other hand, being purely perturbative, it cannot 
generate the mass gap needed to ensure a fast conver- 
gence of the integration over the impact parameter b. 
Because of this problem, up to now all the phenomeno- 
logical applications of the BK equation were based on 
model assumptions regarding the 5-dependence. It is 
always assumed that the 6-dependence factorizes and 
in practice the BK equation is usually solved without 
any trace of b. At the end, the 6-dependence is restored 
via an ansatz with a typically exponential or Gaussian 
profile. An attempt to go beyond this approximation 
has been reported in Ref. [167, 168]. 

— It is very desirable to go beyond the BK equation and 
relax all underlying assumptions outlined above. The 
higher order corrections arc most needed. In particular 
it is important to learn how to include the running of 
as, though in the phenomenological applications the 
running of as is usually implemented. 



^ There was a recent claim of Mueller and Shoshi [176] that 
the s-channel unitarity is in fact violated during the evolution. 
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— The LO BFKL kernel does not have the correct short 
distance limit responsible for the Bjorken scaling vi- 
olation. As a result the BK equation does not natu- 
rally match with the DGLAP equation. Though sev- 
eral approaches for unification of the BK equation 
and DGLAP equations were proposed [52, 162, 165, 
177, 178], the methods are not fully developed. All 
approaches deal only with low x and only with the 
gluon sector. We would like to have a unified evolution 
scheme for both small and large x and with quarks 
included. 



5.2 Phenomenology with the BK equation 

The deep inelastic structure function F2 is related to the 
dipole amplitude N via 



F2ix,Q^) = 



—J / d rx / dzP'^ (Q ; r_L, z)CTdipoic(?'±, a;), 



(53) 



with the dipole cross section given by the integration over 
the impact parameter: 

<Jdipoic{r±,x) = 2 / d'^bN{r^,x;b). (54) 

The physical interpretation of eq. (|53|l is transparent. 
It describes the two stages of DIS [179]. The first stage is 
the decay of a virtual photon into a colorless dipole {qq - 
pair) . The probability of this decay is given by P^ known 
from QED [134,180-182]. The second stage is the interac- 
tion of the dipole with the target (adipoic hi eq. (|53(l ). In 
the large Nc limit a color charge has a well-defined anti- 
charge partner in a color dipole. Eq. (|53|l illustrates the 
fact that in this limit these color dipoles are the relevant 
degrees of freedom in QCD at high energies [134]. 

For the phenomenological applications one may use 
the function N{r±, x; b) or (Tdipoic(?'_L, x) obtained directly 
from the solutions of the BK equation H51|) . With addi- 
tional DGLAP corrections this approach was adopted by 
Gotsman et al. in Ref. [162]. 

Alternatively one can relate N to an unintegrated 
gluon distribution function !F{x,k'^) = f{x,k'^)/k'^. The 
dipole cross section can be expressed via / [183, 184]: 



cr dipoic{r±,x) 



^J^[l-Mkr^)]as{k')f{x,k') 
The inversion of eq. H55() is straightforward 



fix,k^) 



(fbh{k^,x,b), 



h{k ,x,b) 



N, 



A as TT^ 



as TT^ d (In fc2) 



k'^AkN{k^,x,b) 

N{P,x,b). 



(55) 



(56) 



(57) 



Here Ak is the 2-dimensional Laplace operator. The func- 
tion N is related to the Fourier transform of N 



N{k^,x,b) = 



(Pr± 
2 7rr2 



N{r^,x,b). 



(58) 



In fact, N obeys a nonlinear version of the LO BFKL 
evolution equation in momentum space. The function N 
can be interpreted as an unintegrated gluon distribution. 
N and h coincide at large momenta but differ at small 
ones. On one hand, within the dipole picture it is rather 
the function N which gives the probability to find a gluon 
with a given transverse momentum and at a given impact 
parameter. On the other hand, it is the function / (or h) 
which enters the k± (high energy) factorization formula. 
In what follows we will concentrate on the unintegrated 
gluon distribution / only. 

Instead of solving the BK equation H51(l and then in- 
verting the relation (|55|1 one can adopt another strategy 
and reformulate the problem directly in terms of the un- 
integrated gluon density /. This approach was adopted in 
work by Kutak-Kwiecinski [178] and Kutak-Stasto [185]. 
Using relations 154|) and (|55|l one can transform l|51(l into 
an equation for the unintegrated gluon distribution 



f{x,e)^f'^"Hx,k') 



+ 



N^as{P)_^2 f ± 

IT j^ Z 

' X 1.12 



dk^ 
fc'2 



/(f,A:'^)-/(f,A:2) , /(f,fc^) 



\k'^ 
d 



fc2| 
i?2 



|4/i:'4 + fc4|5 



(59) 



^-^^^' -/ - 



1 dz 



fc2 



dk'H'^asik'^)\n{k'^k^) fiz,k'^) 



Here it is written as an integral equation, corresponding 
to the BFKL equation in momentum space supplemented 
by the negative nonlinear term. The input f^^'{x,k'^) is 
given at the scale /cq = IGeV^. This equation was derived 
under the following factorization ansatz: 



N{k^,x,b) = fi{k'^,l,x)S{b) 



(60) 



with normalization conditions on the profile function S{b) 



fd^hSib) = l; fd^hS^ib) 



1 



7ri?2- 



(61) 



The assumption H60|) is crude and corresponds to a situ- 
ation where the projectile size (color dipole) is neglected 
compared to the target size (proton). A simple way to im- 
prove H59|l is to implement NLO corrections in the linear 
term of the equation. It can be done within the unified 
BFKL-DGLAP framework which is presented below. The 
final equation (eq. H64|) below) can be used for phenomeno- 
logical applications. Figs. EOI [2] display the unintegrated 
gluon distributions / obtained in Refs. [162, 185]. 
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Fig. 20. The unintegrated gluon distribution f{x,k'^) as a Fig. 21. The unintegrated gluon distribution f{x,k^) as a 



function of x for different values k^ — 1.5 GeV^ (top) and 
k^ — 30 GeV^ (bottom). Solid lines correspond to the solution 
of the nonlinear equation using GLLM [162] parameterization 
whereas dashed lines (KKS) correspond to the approximate 
solution of I64II [178,185]. For reference we also present linear 
BFKL/DGLAP evolution (KMS) [52] . 



5.3 The saturation scale 



function of k for different values a; = 10 (left) and a; = 10 
(right). The lines are the same as in fig. 1201 . 



In order to quantify the strength of effects that slow down 
the gluon evolution one introduces the saturation scale 
Qs{x). It divides the (x, fc^)-space into regions of dilute 
and dense partonic systems. In the case when fc^ < Ql{x) 
the solution of the BK equation exhibits geometrical scal- 
ing, which means that it depends on one variable only, 
N(r,x) = N{r Qs{x)) or in momentum space N{k'^,x) = 
N{k/Qs{x)). In Fig. [23 we present saturation scales ob- 
tained from H59() in [185] and the corresponding result ob- 
tained from Ref. [162]. Note, however, that the satura- 
tion scale is defined differently in these two approaches. 
In ref. [185] the saturation scale is defined quantitatively 
as a relative difference between the solutions to the linear 
and nonlinear equations, while in ref. [162] it is defined 
by the requirement that N{2/Qs,x) is constant equal to 
1/2 or 1/e. Note that both the KKS [185] and the GLLM 
models predict a saturation scale much bigger than the 
one from the GBW model. 




Fig. 22. Saturation scale from various models. The solid lines 
defines a band of possible saturation scales coming from the 
GLLM model [162]. The dashed line (KKS) is from Ref. [185]. 
The dotted line (GBW) is the Golec-Biernat Wiisthoff model 

[151]. 
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with the scale k^ . The final improved nonlinear equation 
for the unintcgrated gluon density becomes 



Saturation effects are most easily studied in the coordi- 
nate space representation in which it has been difficult 
to include non-leading effects, and the non-leading effects 
have mainly been studied in momentum space, where it is 
hard to include saturation. We now present a short (not 
complete) review of recent theoretical activities which at- 
tempt to go beyond the leading order BK equation. An 
important issue relating to the NLO corrections is energy- 
momentum conservation, which was already addressed in 
section 1^21 and will be further discussed in more detail in 
section IH31 



5.4.1 Beyond leading order 



T^ L z Jki k'^ 



asjk^ 

27r 



|/fc'2-/c2| 



dzPgg{z 



' e 



|4/c'4-|-fc4|5 






_ (l_fc^^) ^ / ^ 



dk^ j i?2 

r-OO 77/2 



r^-»""°(S)^"-'"* 



(64) 



with the input distribution /(°)(a;, A;^). 



The BFKL kernel is known at next-to-leading order. Nev- 
ertheless, a nonlinear equation at NLO has not been de- 
rived yet. I. Balitsky and A. Belitsky [186] have been able 
to compute a single NLO contribution which has maximal 
nonlinearity, namely the N^ term: 



dN 

—— = asKer® {N 
ay 



N-N)~a'^^Ker(g,N-N-N. (62) 



The new kernel Ker can be found in Ref [186]. Triantafyl- 
lopoulos [187] has considered NLO BFKL in the presence 
of a saturation boundary. The results show a decrease 
in the saturation scale growth as a function of rapid- 
ity towards the value A ~ 0.3 observed experimentally 
(GBW [151] and GLLM [162] models). 

Another approach [177, 178] to partially include the 
NLO corrections into the BK equation is to implement 
in the hnear term of eq. (EHl the unified BFKL-DGLAP 
framework developed in [52]. In this scheme the BFKL 
kernel also gets modified by the consistency constraint [17, 
24, 188] 



fc'2 < k^z. 



(63) 



The origin of this constraint is the requirement that the 
virtuality of the exchanged gluon is dominated by its 
transverse momentum [/c'^j ~ k'^ (see also section IT^ . 
The constraint (|63|l resums a large part of the subleading 
corrections in Inl/a;, and it is also connected to the con- 
servation of the negative lightcone component p_ = E—pL 
{cf. section lF31l . Additionally, the non-singular part of the 
leading order DGLAP splitting function, which infiuences 
the normalization of the unintcgrated gluon distribution, 
is included into the evolution and as is assumed to run 



5.4.2 JIMWLK 

The Nc corrections can be accounted for through the 
JIMWLK functional equation [189-193], which is equiva- 
lent to Balitsky's original infinite chain of equations [147]. 
Introducing iV as a target expectation value of a certain 
operator (product of two Wilson hues) , A'' = {W) target, 
the first couple of equations of the Balitsky chain are 



d{W] 
dy 



diWW] 



a^Ker ® {{W) - {WW)), 



(65) 



dy 



asKer ® {{WW) - {WWW)). (66) 



The large Nc limit and the absence of the target corre- 
lations used by Kovchegov [144] is equivalent to a mean 
field approximation which allows to express a correlator 
of a product as a product of correlators: 



WW) = {W){W) = NN- Nr 



oo. 



Thus the first equation of the Balitsky chain closes to the 
BK equation. 

Rummukaincn and Weigert [166] have produced a first 
numerical solution of the JIMWLK equation. They do not 
find any qualitative deviation from solutions of the BK 
equation. The Nc corrections were found to be at a level 
of a few percents. 

Bartels, Lipatov, and Vacca [149] have considered Nc 
corrections to the triple Pomeron vertex: 



dN 
dy 



Us Ker (g) {N - N N 



1 



where the function n has to satisfy a separate equation. 
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5.4.3 Target correlations 

For proton and realistic (not very dense) nucleus targets 
a systematic approach towards inclusion of target corre- 
lations has been developed by Levin and Lublinsky [194]. 
Target correlations can be introduced via a certain lin- 
ear functional differential equation. In general, this linear 
functional equation cannot be reformulated as a non-linear 
equation. However, in a particular case when all n-dipole 
correlations can be accounted for by a single correlation 
parameter, the equation can be brought to a modified ver- 
sion of the BK equation: 

dN 

= asKer ® {N - kNN); (67) 

ay 

In Eq. H67|l k > 1 is the correlation parameter to be found 
from a model for the target. 



5.4.4 Pomeron loops 

Pomeron loops are the first steps towards restoration of 
the i-channel unitarity. lancu and Mueller [195] have con- 
sidered rare fluctuations which were interpreted by Kozlov 
and Levin [196] as pomeron loop contributions. Unfortu- 
nately, it looks as if contributions from the pomeron loops 
are difhcult to incorporate in a framework of a single equa- 
tion. They are known to modify the asymptotic behavior 
of the amplitude N in the deep saturation limit, where 
they give the following asymptotic behavior: 



N{Y) 
N{Y) 



-c{Y-Yo)\ 



Y 



-c(y-Yo)V2. 



2 a, BKE 



Y —> oo Pom Loops 



Recently there has been a lot of activity in attempting 
to consistently include Pomeron loops into high energy 
evolutions [197-202]. 



5.4.5 Local multi-pomeron exchange 

It is claimed that the BK equation sums all possible con- 
tributions which are not suppressed either by ag or Nc- 
For example, the cubic term in Eq. H62|l appears at next- 
to-leading as order only. In particular it is implied that all 
multi-pomeron exchanges and multi-pomeron vertices are 
either absorbed by the triple pomeron vertex of the BK 
equation or suppressed. Levin and Lublinsky [194] have 
argued that this might be not true. They argue that in 
addition to a possibility for a pomeron to split into two, 
there exists a process of multi-pomeron exchange, which is 
local in rapidity. After these contributions were resummed 
in the eikonal approximation, a new modification of the 
BK equation was proposed: 



dN 
dy 



—— = {\ ~ N)as Ker (X) {N - N N). 



(68) 



5.5 Energy conservation aspects 

Main author G. Gustafson 



5.5.1 Rapidity veto 

It is well known [16] that a major fraction of the higher 
order corrections to (not only) BFKL is related to energy 
conservation. The large effect of energy-momentum con- 
servation is also clearly demonstrated by the numerical 
analyses by Andersen-Stirling [21] and Orr-Stirling [18]. 
Conservation of energy and momentum implies the con- 
servation of both the positive and the negative lightcone 
components, p± = E ±pL. Although most analyses have 
concentrated on the conservation oi p+, as being more im- 
portant, we will see below that also conservation of p_ has 
a very significant effect. In LLA the steps in ln(l/a;) are as- 
sumed to be large, and the necessary recoils due to energy 
conservation are neglected. The main effect of conserva- 
tion of the positive lightcone component p+ = E + p^, is 
that small steps in ln(l/a::) with corresponding large recoils 
are suppressed. One way to take this into account is to in- 
troduce a veto, not allowing steps in ln(l/x) smaller than 
a cut rj. (This is called a rapidity veto also if the evolution 
variable is defined as y = ln(l/a;) and not the true rapid- 
ity.) The effect of such a veto is studied in refs. [75,78,203], 
and at high energies it has a similar effect as the higher 
order corrections, reducing the growth at small x. 

A recent study of the BK equation in the presence of 
a rapidity veto is presented by Chachamis, Lublinsky and 
Sabio Vera [204]. The application of this method to the 
BK equation makes it non-local in rapidity: 



dNjy) 
dy 



= ttsKer (g) {N{y-r]) ~ N{y-ri) N{y-T])). 



The veto somewhat delays saturation in accordance with 
the expectations associated with the next-to-leading order 
corrections. If the veto is put on top of the BK equation 
with running as then the effect of NLO corrections is sig- 
nificantly reduced. This observation gives support to the 
phenomenological studies of Refs. [162, 178]. 

An similar approach to this problem is presented by 
Gotsman, Levin, Maor, and Naftali [205]. The effects of 
the cut in ln(l/a;) is taken into account in a modified BK 
equation: 



dN{r,Y;b) Cpas f d^r'r 



dY 



72^/^2 



(r " r')" 



_d_ 
W 



2N { r', Y;h -(r - r') ] - iV (r, Y; b) 
-N{r',Y;h~^{r-r')]N(r~v',Y;h-^r' 



(69) 



The derivative under the integral is related to a cut in 
ln(l/a;) ex lnp+. The modification of the pole at 7 = 1, 
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Fig. 23. A quark-antiquark dipole in transverse coordinate 
space is split into successively more dipoles via gluon emission. 



which is related to conservation of the negative lightcone 
component p_ — E ~ p^ (or the inverse k± ordering) and 
the consistency constraint, is not included. The motivation 
for this is that this effect is not important once the dipole 
density has reached saturation, that is for x so small that 

Q?(^) > Q^- 

We note, however, that the non-leading effects can sig- 
nificantly reduce the value of Q^ (x) , and thus delay the on- 
set of saturation, as discussed in e.g. refs. [187]. An essen- 
tial result of the analysis discussed in the next subsection 
is that also the conservation of p- has an important ef- 
fect, and contributes significantly to pushing the a;- values, 
where saturation becomes essential, to smaller values. We 
note also that an estimate of the relative importance of 
saturation and non-leading effects for the reduced growth 
rate is very important for reliable extrapolations to higher 
energies at LHC and high energy cosmic rays. 



5.5.2 Full energy-momentum conservation 

A different approach to energy-momentum conservation 
is presented in ref. [206]. As discussed above non- leading 
effects are most easily studied in momentum space, while 
unitarity or saturation effects are easier analyzed when 
formulated in transverse coordinate space. In ref. [206] 
similarities between the Linked Dipole Chain model 
(LDC) [203, 207] in momentum space and the Mueller 
dipoles in transverse coordinate space [134, 135, 154] are 
used to derive a scheme for implementing energy momen- 
tum conservation in Mueller's dipole formalism. It is con- 
jectured that only those gluon emissions, which satisfy 
energy-momentum conservation, can correspond to real 
final state gluons, and that keeping only these (with a 
corresponding modification of the Sudakov form factor) 
will not only give a better description of the final states, 
but also account for essential parts of the NLO corrections 
to the BFKL equation. The approach is based on the ob- 
servation that the emission of a dipole with a very small 
transverse size, r, corresponds to having two very well lo- 
calized gluons, and such gluons must have large transverse 
momenta of the order p±^ '^ 1/r. By in this way assign- 
ing a transverse momentum to each emitted gluon, and 
also taking into account the recoils of the emitting glu- 
ons, it is possible to make sure that each dipole splitting 
is kinematically allowed. 



Formalism In the process 7* -^ QQ — + QgQ — + QggQ ~> 
. . ., a virtual photon is split into a QQ color dipole, which 




Fig. 24. A dipole cascade, where a chain of smaller and smaller 
dipoles is followed by a set of dipoles with increasing sizes. 
This is interpreted as one fc_L-ordered cascade from the left 
and one from the right, up to a central hard subcollision, which 
is represented by the dipole with minimum size and therefore 
maximum k±. 



is first split into two dipoles by the emission of a gluon, 
then into three dipoles by a second gluon, etc. The pro- 
cess is illustrated in transverse coordinate space in fig. 1231 
The probability for such a dipole splitting is given by the 
expression [134,135,154] (for notation see fig. 123(1 



dP a 2 

dy 27r i 

S — exp 
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Y2- 
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'02 ' 12 



(70) 



Here S denotes a Sudakov form factor. We note that the 
integral over d'^Y2 in the exponent diverges for small values 
of ro2 and ri2. Therefore Mueller introduced a cutoff p, 
such that the integration region satisfies ro2 > p and ri2 > 
p. A small cutoff value p will here imply that we get very 
many dipoles with small r-values. 

If a dipole size, r, is small, it means that the gluons are 
well localized, which must imply that transverse momenta 
are correspondingly large. This implies that not only the 
new gluon gets a large fc^ ~ 1/r, also the original gluon, 
which is close in coordinate space, gets a corresponding 
recoil. Let us study the example in fig. |^ For the emis- 
sions of the gluons marked 2, 3, and 4 the dipole sizes 
become smaller and smaller, a ^ 6 ^ c ^ d, in each step 
of the evolution. The corresponding k± therefore become 
larger and larger in each step. After the minimum dipole, 
with size d, the subsequent emissions, 5, and 6, give again 
larger dipoles with correspondingly lower fcj^ values. The 
probability for this chain is proportional to 



d^r2 a^ d^rsb^ d^r4C^ d'^r^d'^ d^r^e^ 



152 



62, 



c2d2 



PP 



(71) 



For the first emissions, 2 and 3, in this expression we 
recognize the product of factors d'^Vi/rf ~ Yld'^ki/k^, 
just as is expected from a "DGLAP evolution" of a chain 
with monotonically increasing k± . Emission number 4 cor- 
responds to the minimum dipole size, d, and we here 
note that the factors of d cancel in eq. H71|) . We there- 
fore get the weight d'^r4 ^ d^]i^i:,x/k'^^^, which corre- 
sponds to a hard gluon-gluon collision. When the dipole 
sizes get larger again, this gives factors corresponding to 
a "DGLAP chain" from the other end of the chain, up to 
the central hard subcollision. 

It is also easy to see that for a chain with increasing 
dipole sizes up to a maximum value, Tmax, which thus 
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Fig. 25. A symbolic picture of a 7*7* collision in rapidity- 
rx-space. The two dipole chains interact and recouple with 
probability / given by eq. 17211 . 



corresponds to a mininium transverse momentum, k±„ 



we get the weight d r^^y^/r\ 



,4 

max 



d kmin- Therefore there 



is no singularity for the minimum A:^-value. This result 
agrees exactly with the result in the Linked Dipole Chain 
model, LDC [203,207], which is a reformulation of the 
CCFM model [24,25], interpolating between DGLAP and 
BFKL for non-fcj^-ordered chains. 

To study 7*7* scattering we imagine that the two vir- 
tual photons split up into quark-antiquark pairs, which 
develop into dipole cascades as schematically illustrated 
in fig. 1211 When the two central dipoles collide and inter- 
act via gluon exchange, it implies a recoupling of the color 
charges, as indicated by the arrow, and the probability for 
this is given by the expression [208] 



/^f^ln 



|ri -rsl • |r2 -r4| 



ri -r4 • r2 -ra 



(72) 



As the dipole cascades from the two virtual photons 
branch out, it is also possible to have multiple collisions, 
when more than one pair of dipoles from the left and the 
right moving cascades are interacting. The total cross sec- 
tion is then given by 



(fb{l 



-Ek 



(73) 



where b denotes the impact parameter. 

With a small cutoff p [r > p) we get, as mentioned 
above, very many small dipoles. If these are interpreted 
as real emissions, it would imply a violation of energy- 
momentum conservation. The emission of these small 
dipoles must be compensated by virtual emissions. Thus 
the result in eq. H73|) will describe the inclusive cross sec- 
tion, but the many dipoles produced in all the branching 
chains will not correspond to the production of exclusive 
final states. 

The main feature of the LDC model is the observa- 
tion that both the total cross section and the final state 
structures are determined by chains consisting of a subset 
of the gluons appearing in the final state. These gluons 
were called "primary gluons" in ref. [203] and later called 
"backbone gluons" in ref. [209] . Remaining real final state 
gluons can be treated as final state radiation from the 
primary gluons. Such final state emissions do not mod- 
ify the total cross sections, and give only small recoils to 



the parent emitters. The primary gluons have to satisfy 
energy-momentum conservation, and are ordered in both 
positive and negative light-cone momentum components, 
p+ and P-. We saw above that in Mueller's cascade the 
emission probabilities for gluons, which satisfy the condi- 
tions for primary gluons in LDC, have exactly the same 
weight, when the transverse momenta are identified with 
the inverse dipole size, 2/r. This inspires the conjecture 
that with this identification an appropriate subset of the 
emissions in Mueller's cascade can correspond to the pri- 
mary gluons in the momentum space cascade, meaning 
that they determine the cross sections while the other 
emissions can be regarded as either virtual fiuctuations 
or final state radiation. 

A necessary condition for this subset of gluons is that 
energy and momentum is conserved. Therefore we ex- 
pect that keeping only emissions which satisfy energy- 
momentum conservation can correspond to real emissions, 
and keeping only these emissions (with a corresponding 
modification of the Sudakov form factor) will not only ac- 
count for important NLO effects, but also give a closer 
correspondence between the generated dipole chains and 
the observable final states. 

A very important consequence of energy-momentum 
conservation is also that it implies a dynamical cutoff, 
p{Ay), which is large for small steps in rapidity. Ay, but 
gets smaller for larger Ay. (Alternatively it could be de- 
scribed as a cutoff for Ay which depends on r. Note that 
in this formalism y is the true rapidity and not log(l/a;).) 
Conserving also the negative light-cone momentum, p_, 
implies that in a similar way we may also get a maximum 
value for r in each emission. 

The net result of conservation of both p+ and p- 
is that the number of dipoles grows much more slowly 
with energy. Besides its physical effects, this also sim- 
plifies the implementation in a MC program. It is here 
straight forward to calculate cross sections and to study 
saturation effects, by comparing the unitarized expression 
J (fb{l — e~^ f^^ ) in eq. ifT^ with J d'^bj^ ftj representing 
single IP exchange. (The large numerical complications in 
MCs without energy conservation, discussed in ref. [208], 
are not present.) 



Results Below we show some results obtained with a fixed 
coupling a — 0.2. 

Dipole- dipole scattering. The cross section for scatter- 
ing of two dipoles with sizes ri and r2 is shown in fig. 1261 
With a fixed coupling the scaled cross section, o-/r|, de- 
pends only on the ratio ri/r2. We can imagine a tar- 
get with size r2 ^ ^/M, and a varying projectile size 
fi ~ Xj \fQ^. The results show that the cross section 
grows faster with the total rapidity range, y ~ Ins, for 
smaller r\ (larger Q^), in a way qualitatively similar to 
the behavior of the proton structure function. 

The effect of energy conservation is demonstrated in 
fig.|57|by the results obtained for the case r\ =r2, with a 
constant cutoff, p — 0.02 r^. Comparing with fig.EHlwe see 
that energy conservation has a very strong effect, reducing 
(J by almost an order of magnitude for 1^ ~ 13. 
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Fig. 26. The scaled unitarized dipole-dipole cross section, 
cr/r2, as a function of Y for different initial conditions. 
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Fig. 27. The scaled unitarized (full line) and one-pomeron 
(dashed line) dipole-dipole cross sections calculated without 
energy conservation. 



In fig. 1271 we also see that without energy-momentuiii 
conservation the effect of multiple IP exchange (satura- 
tion) is about a factor 2 for ri = r2 and Y = 13. 
The much smaller cross section obtained with energy- 
momentum conservation implies that the saturation effect 
is much less important, being only « 20% for the same pa- 
rameter values. 

Dipole-nucleus and dipole-proton collisions. Dipole- 
nucleus collisions have been studied using a toy model 
nucleus, with a Gaussian distribution in dipole size r and 
impact parameter b. The dipole density is given by 



dN ■ 



B-d^re" 



"^'/'■o.d2be-''Vfc^ 



(74) 



The widths of the distributions are taken to be rg = 1 fm 
and bo = A^^^ ■ 1 fm (where A is the mass number of the 
nucleus), and the normalization constant B is adjusted so 
that the transverse energy is given by A-1 GeV. 

The results for A — 200 and projectile sizes rproj =0.1 
and 1 GeV~^ are shown in figure 1^ Results are presented 
both for single pomeron exchange and including unitariza- 
tion. The effect of unitarization grows with nuclear size 
and with the size of the projectile. For a small projectile 




Fig. 28. The dipole-nucleus cross section for rproj =0.1 and 
1 GeV~ and A = 200. The unitarized result is shown by the 
solid lines, and the one-pomeron contribution by the dashed 
lines. 



of size 0.1 GeV^^ we can see the effect of color trans- 
parency, as the cross sections for the unitarized and the 
one pomeron calculations are almost identical. For a larger 
projectile we do see a clear effect from unitarization, but 
even for rp^oj = 1 GeV~^ and a nucleus with A = 200 this 
effect is only about 20% in the rapidity interval 10— 14. For 
smaller nuclei the effect will be correspondingly smaller. 

When the same toy model is applied to deep inelastic 
ep scattering (with A — 1 and simply identifying Q^ with 
4/rpj.Qj), we want to emphasize that we here only want to 
study the qualitative behavior. A quantitative comparison 
with HERA data has to wait for an improvement of the 
crude toy model for the proton target (dipole correlations 
may be important), and one should then also take into 
account the detailed effects of the photon wavefunction. 

The resulting dipole-nucleon cross section is shown in 
figure 123 for two different projectile sizes, corresponding 
to g2 = 4 GeV^ and Q^ = 400 GeV^ The result for single 
pomeron exchange, i.e. without unitarization corrections, 
is shown by the dashed lines, and we see that the effect 
from unitarization is quite small. 

In figure 123 we also see that the logarithmic slope 
Acff = d{\og a)/d{log 1/x) is increasing with increasing Q^. 
The effective slope, Aeff , is not a constant for fixed Q^, but 
depends on both Q^ and x, when unitarization and/or 
energy conservation is taken into account. For the com- 
parison with experimental data figure 1201 shows Acff deter- 
mined in the x-interval used in the analysis by HI [210], 
which varies from a; « 2 x 10~^ for Q^ — 1.5 GeV^ to 
X w 3 X 10-2 for Q2 = 90 GeV^. We note that the result 
of our crude model is not far from the experimental data, 
although the dependence on Q^ is somewhat weaker in the 
model calculations. As in figure H^ we see that the effect 
of unitarization is small, and, as expected, it gets further 
reduced for larger Q^- values. 

Thus we find that the result of the simple model is 
surprisingly close to experimental data from HERA. The 
effect of energy conservation is a suppression for small x- 
values and small Q^, which is qualitatively similar to the 
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Fig. 29. The scaled dipole-p cross section as a function of 
logl/a;, for Q^ = 4 GeV^ and Q^ = 400 GeV^ The unitarized 
results are shown by the solid lines while the dashed lines show 
the one-pomeron results. 
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Fig. 30. The effective slope measured at different Q^ com- 
pared to data from HERA. The full line is our model including 
unitarization, while the dashed line is without. Filled circles 
are data from ZEUS [211], filled [31] and open [212] squares 
are data from HI. 



effect expected from unitarization. This suppression is so 
strong that the effect from adding unitarization is only a 
very small correction, visible for small (^--values. 

If we compare these results with those of ref . [205] , we 
find a significantly larger effect from energy-momentum 
conservation. One reason appears to be the inclusion of 
p_ -conservation. This is related to the consistency con- 
straint in eq. (|6!-{|l . which orders the emissions in the nega- 
tive lightcone momentum. In the formalism discussed here 
this is found to have a noticeable effect. Thus we find that 
including only conservation of p+, and not of p_, increases 
the cross section by a factor 2 (3) for dipole-proton col- 
lisions at Q^ = 4 (400)GeV^. Consequently we conclude 
that full energy-momentum conservation is very essential 
for the result and for the relative importance of saturation 
and NLO effects. 



5.6 Outlook 

It is essential for the future phcnomenological studies to 
eliminate the model dependent treatments of the impact 
parameter. Though the BK equation has been solved nu- 
merically with the full 5-dependence traced [167, 168], 
these results are not yet suitable for phenomenological ap- 
plications. 

A further study of the relation between the dipole pic- 
ture vs. traditional diagrammatics based on the s-channel 
unitarity is needed. In particular, it is not clear if the 
dipole picture survives at NLO. In general there is a quest 
for a simple effective Reggeon field theory in QCD. 

The large effect of full energy-momentum conservation 
make further studies of the relative importance of NLO 
effects and saturation important. 

NLO effects and saturation both contribute to a reduc- 
tion of the parton distributions for small x. An improved 
understanding of these effects, including the relation be- 
tween them, is very important for extrapolations to higher 
energies at LHC or high energy cosmic ray events. 

The discussions presented above concentrate on total 
or inclusive cross sections. More work is also needed to 
calculate the properties of the resulting final states. 



6 Multiple interactions, saturation and 
rapidity gaps 

6.1 AGK cutting rules 

Main author J. Bartels 



6.1.1 Introduction 

Multiple parton interactions play an important role both 
in electron proton scattering at HERA and in high energy 
proton proton collisions at the LHC. At HERA, the linear 
QCD evolution equations provide, for not too small Q^, a 
good description of the F2 data (and of the total 7*7? cross 

section, cr^ot)- This description corresponds to the emis- 
sion of partons from a single chain (Fig. I31h .l. However, 
at low Q^ where the transition to nonperturbative strong 
interaction physics starts, this simple picture has to be 
supplemented with corrections. First, there exists a class 
of models [151, 153, 213] which successfully describe this 
transition region; these models are based upon the idea of 
parton saturation: they assume the existence of multiple 
parton chains (Fig. ISlb ) which interact with each other, 
and they naturally explain the observed scaling behavior, 
^2(Q',2:) « F2{Q^/Ql{x)) ^ii\iQl{x) = Ql(l/x)^. Next, 
in the photoproduction region, Q^ sa 0, direct evidence 
for the presence of multiple interactions also comes from 
the analysis of final states [214]. A further strong hint at 
the presence of multi-chain configurations comes from the 
observation of a large fraction of diffractive final states in 
deep inelastic scattering at HERA. In the final states anal- 
ysis of the linear QCD evolution equations, it is expected 
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Fig. 31. Contributions to the total cross section a'^^f: (a) the 
single chain representing the linear QCD evolution equations; 
(b) gluon production from two different gluon chains. 



WW 





Fig. 32. Hard diffractive final states. (a) dijet production; (b) 
the diffractive cross section as s-channel discontinuity of a two- 
ladder diagram. 



that the produced partons are not likely to come with 
large rapidity intervals between them. In the momentum- 
ordered single chain picture (Fig. Et^), therefore, diffrac- 
tive final states should be part of the initial conditions 
(inside the lower blob in Fig. 13 Ih ,). i.e. they should lie 
below the scale Ql which separates the parton descrip- 
tion from the nonperturbative strong interactions. This 
assignment of diffractive final states, however, cannot be 
complete. First, data have shown that the Pomeron which 
generates the rapidity gap in DIS diffraction is harder than 
in hadron - hadron scattering; furthermore, there are spe- 
cific diffractive final states with momentum scales larger 
than Qo, e.g. vector mesons built from heavy quarks and 
diffractive dijets (illustrated in Fig. I32|l : the presence of 
such final states naturally requires corrections to the sin- 
gle chain picture (Fig. I82b'l. From a i-channel point of 
view, both Fig. 131b and Fig. 132b belong to the same class 
of corrections, characterized by four gluon states in the 
i-channel. 

In proton-proton collisions corrections due to multi- 
ple interactions should be important in those kinematic 
regions where parton densities for small momentum frac- 
tions and for not too large momentum scales are being 
probed, e.g. jet production near the forward direction. 
Another place could be the production of multijet final 
states (Fig. \^ : multiple jets may come from different 
parton chains, and these contributions may very well af- 
fect the background to new physics beyond the standard 




Fig. 33. Jet production in pp collisions from two different 
parton chains 



model. Moreover, the modeling of multijet configurations 
will be necessary for understanding the underlying event 
structure in pp collisions (see [215] and references therein). 

From the point of view of coUinear factorization, multi- 
ple interactions with momentum ordered parton chains are 
higher- twist effects, i.e they are suppressed by powers of 
the hard momentum scale. At small x, however, this sup- 
pression is compensated by powers of the large logarithms, 
Inl/x: multiple interactions, therefore, are mainly part of 
small-x physics. In this kinematic region the Abramovsky- 
Gribov-Kanchelli (AGK) [216] rules can be applied to the 
analysis of multi-gluon chains, and it is the aim of this ar- 
ticle to present a brief overview about the current status 
of the AGK rules in pQCD. 

As we will discuss below, in the analysis of multiple 
parton chains the couplings of n gluons to the proton play 
an essential role. Regge factorization suggests that these 
couplings should be universal, i.e. the couplings in 7*p 
collisions at HERA are the same as those in pp scatter- 
ing at the LHC. Therefore, a thorough analysis of the role 
of multiple interactions in deep inelastic electron-proton 
scattering at HERA should be useful for a solid under- 
standing of the structure of events at the LHC. 



6.1.2 Basics of the AGK cutting rules 

The original AGK paper [216], which was written before 
the advent of QCD, addresses the question how, in the 
optical theorem. 






(75) 



the presence of multi-Pomeron exchanges (Fig. I34|l in 
the total hadron-hadron cross section leads to observ- 
able effects in the final states (rhs of ea. (j75|l ). Based 
upon a few model-independent assumptions on the cou- 
plings of multi-Pomeron exchanges to the proton, the 
authors derived simple 'cutting rules': different contri- 
butions to the imaginary part belong to different cuts 
across the multi-Pomeron diagrams, and each cut has its 
own, quite distinct, final state characteristics. As a re- 
sult, the authors found counting rules for final states with 
different particle multiplicities, and they proved cancella- 
tions among rescattering corrections to single-particle and 
double-particle inclusive cross sections. 
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Fig. 35. Decomposition of the coupling of four gluons to a 
virtual photon. In the last two terms on the rhs it is understood 
that we have to sum over different pairings of gluons at the 
lower end. 



Fig. 34. s-cut through a multi-Pomeron exchange: the zig-zag 
lines stand for nonperturbative Pomerons. 



In the QCD description of hard (or semihard) final 
states a close analogy appears between (color singlet) 
gluon ladders and the nonperturbative Pomeron: multiple 
parton chains (for example, the two chains in Fig. 131b ) 
can be viewed as cuts through two perturbative BFKL 
Pomerons. In the same way as in the original AGK pa- 
per, the question arises how different cuts through a QCD 
multi-ladder diagram can be related to each other. In 
the following we briefly describe how AGK cutting rules 
can be derived in pQCD [217,218]. Subsequently we will 
present a few new results which come out from pQCD cal- 
culations, going beyond the original AGK rules, followed 
by some numerical estimates of the effects which can be 
expected. 

One of the few assumptions made in the original AGK 
paper states that the coupling of the Pomerons to the ex- 
ternal particle are (i) symmetric under the exchange of 
the Pomerons (Bose symmetry), and (ii) that they remain 
unchanged if some of the Pomerons are being cut. These 
properties also hold in pQCD, but they have to be refor- 
mulated: (i') the couphng of (reggeized) gluons to exter- 
nal particles is symmetric under the exchange of reggeized 
gluons, and (ii') it remains unchanged if we introduce cut- 
ting lines between the gluons. In QCD, however, the color 
degree of freedom also allows for another possibility: in- 
side the n-giuon state (with total color zero), a subsys- 
tem of two gluons can form an antisymmetric color octet 
state: in this case the two gluons form a bound state of 
a reggeized gluon (bootstrap property). For the case of 
7*7* scattering, explicit calculations [219] have shown 
that the coupling of n gluons to virtual photons can be 
written as a sum of several pieces: the fully symmetric 
('irreducible') one which satisfies (i') and (ii'), and other 
pieces which, by using the bootstrap property, can be re- 
duced to symmetric couplings of a smaller number of glu- 
ons ('cut reggeons'). This decomposition is illustrated in 
Fig. 1211 Since the bootstrap property is related to the 
reggeization of the gluon and, therefore, is expected to be 
valid to all orders of perturbation theory, also these prop- 
erties of the couplings of multi-gluon states to external 
particles should be of general validity. In this short review 
we will mainly concentrate on the symmetric couplings. 

As an illustrative example, we consider the coupling of 
four gluons to a proton. The simplest model of a symmet- 




Fig. 36. The symmetric coupling of four gluons to an external 
particle. The lines inside the blob denote the color connection, 
e.g. the first term has the color structure SaiaiSa-^ai- 
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Fig. 37. Different cutting lines in the four-gluon exchange. 



ric coupling is a sum of three pieces, each of which contains 
only the simplest color structure: The best-known cutting 
rule for the four gluon exchange which follows [217,218] 
from this symmetry requirement is the ratio between the 
three different pairings of lines given in Fig.|^ Each term, 
on the partonic level, corresponds to a certain multiplicity 
structure of the final state: a rapidity gap ('zero multiplic- 
ity'), double multiplicity, and single multiplicity. Simple 
combinatorics then leads to the ratio [216] 



1 : 2 : -4. 



(76) 



for the two-ladder contribution to the cross section. In 
order to be able to generalize and to sum over an arbitrary 
number of gluon chains, it is convenient to use an eikonal 
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Fig. 38. AGK cancellations in the one-jet inclusive cross sec- 
tion. 




Fig. 40. Decomposition into two rapidity intervals: the upper 
(left) interval has double multiplicity, the lower (right) one 
corresponds to a rapidity gap. 



to include also the soft Pomeron: all one needs to assume 
is that the couplings of soft Pomerons and reggeized glu- 
ons are symmetric under interchanges, and they are not 
altered if cutting lines are introduced. 



Fig. 39. (a) Nonvanishing rescattering corrections in the one- 
jet inclusive cross section; (b) a new vertex: g + 2g —* jet. 



ansatz: 
N^Jki,ai;.. . ;fc2n, a2„; tj) 



x/W 



=f==f Yl (I)'^{kl,k2;uji2)6a,a2- ■■■ 
c / \ Pairings 



{^2n-l^ k2n\ ^2n-l,2n)^a2„-ia2 



(77) 



Inserting this ansatz into the hadron-hadron scattering 
amplitude, using the large- iVc approximation, and switch- 
ing to the impact parameter representation, one obtains, 
for the contribution of k cut gluon ladders, the well-known 
formula: 



where 



ImAk = 4s / Sbe^^^P{s, b) 



P{s,b) 






k\ 



(78) 



(79) 



and J? stands for the (cut) two-gluon ladder. 

Another result [218] which follows from the symmetry 
properties of the n gluon-particle coupling is the cancel- 
lation of rescattering effects in single and double inclu- 
sive cross sections. In analogy with the AGK results on 
the rescattering of soft Pomerons, it can be shown that 
the sum over multi-chain contributions and rescattering 
corrections cancels (Fig. I38|l , leaving only the single-chain 
contribution (in agreement with the factorization obtained 
in the coUinear analysis). This statement, however, holds 
only for rescattering between the two projectiles: it does 
not affect the multiple exchanges between the tagged jet 
and the projectile (Fig. 1^ which require a separate dis- 
cussion (see below). All these results can be generalized 



6.1.3 New results 

Explicit calculations in QCD lead to further results on 
multiple interactions. First, in the four gluon exchange 
there are other configurations than those shown in Fig. 1371 
one example is depicted in Fig. QHI Here the pairing of 
gluon chains switches from (14) (23) in the upper part (= 
left rapidity interval) to (12) (34) in the lower part (= right 
rapidity interval). One can show that the ratio 1:2: —4 
holds for each rapidity interval. In [218] this has been 
generalized to an arbitrary number of exchanged gluon 
lines. 

Another remark applies to the applicability of the cut- 
ting rules to rescattering corrections in the single jet in- 
clusive cross section (Fig. I39|) . Below the jet vertex we, 
again, have an exchange of four gluon lines, similar to the 
diagram in the middle of Fig. |^ As to the cutting rules, 
however, there is an important difference between the two 
situations. In Fig. 133 the blob above the four gluons is 
totally inclusive, i.e. it contains an unrestricted sum over 
s-channel intermediate states, whereas in Fig. 1391 the part 
above the four gluon state is semi-inclusive, i.e. it con- 
tains the tagged jet. This 'semi-inclusive' nature destroys 
the symmetry above the four gluon states, and the cutting 
rules have to be modified [220,221]. In particular, eqs.jZJj) 
- H78|) are not applicable to the rescattering corrections 
between the jet and projectile. A further investigation of 
these questions is in progress [222]. 

Finally a few comments on reggeization and cut 
reggeons. Clearly there are more complicated configura- 
tions than those which we have discussed so far; an ex- 
ample appears in 7*p scattering (deep inelastic electron 
proton scattering). In contrast to pp scattering, the cou- 
pling of multi-gluon chains to the virtual photon can be 
computed in pQCD, and the LO results, for the case of 
n = A gluons, are illustrated in Fig. ^2 It turns out that 
we have two alternative possibilities: in the completely in- 
clusive case (total cross section), it is convenient to chose 
Fig. 141b . i.e. the sum of all contributions can be decom- 
posed into two sets of diagrams. In the first set, at the top 
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Fig. 41. Four-gluon contributions to 7*p proton scattering: 
two equivalent ways of summing over all contributions, (a) 
the decomposition of Fig. 1351 with the pQCD triple Pomeron 
vertex, (b) an alternative way of summation which explicitly 
shows the coupling of two Pomerons to the photon vertex and 
which leads to a new vertex Z. 



of the diagram two gluons couple to the quark- antiquark 
pair, and the subsequent transition to the four-gluon state 
goes via the pQCD triple Pomeron vertex. This vertex, as 
a function of the 4 gluons below, has the symmetry proper- 
ties described above. As a result, we can apply the cutting 
rules to the four gluon state, as discussed before. However, 
there is also the second term in Fig. 141b. which consists of 
a two gluon state only: this is the reggeizing contribution 
we have mentioned before. As indicated in the figure, the 
splitting of the reggized gluons at the bottom amounts to 
a change in the (nonperturbative) coupling. We want to 
stress that, because of the inclusive nature of this set of 
diagrams, the triple Pomeron vertex V in Fig. 141b .. similar 
to the BFKL kernel, contains both real and virtual con- 
tributions. For this reason, the decomposition in Fig. 141b 
is applicable to inclusive cross sections, and it is not con- 
venient for investigating specific final states such as, for 
example, diffractive final states with a fixed number of 
quarks and gluons in the final state. 

There exists an alternative way of summing all con- 
tributions (Fig. 141b ) which is completely equivalent to 
Fig. 141b but allows to keep track of diffractive qq, 
qqg,- ■ ■ final states: this form is illustrated in Fig. 141b. One 
recognizes the 'elastic intermediate state' which was not 
visible in Fig. 141b . and the new triple Pomeron vertex Z 
which contains only real gluon production. This vertex Z, 
as discussed in [223] is no longer symmetric under permu- 
tations of the gluons at the lower end; consequently, we 
cannot apply the AGK cutting rules to the four gluon 
states below. These findings for multiple scattering ef- 
fects in DIS imply, strictly speaking, that cross sections 
for diffractive qq or qqg states cannot directly be inserted 
into the counting rules (|76() . 

Also pp scattering will contain corrections due to mul- 
tiple interactions which are more complex. There are, for 
example, graphs which contain the 2 — > 4 gluon vertex V^ 





Fig. 42. A correction in which the number of lines changes. 
The black vertex denotes the 2^4 gluon vertex. 



leading to a change of the number of gluon lines fFig.l42|l. 
Since this 2^4 gluon vertex, as a function of the four 
gluons below the vertex, satisfies the symmetry require- 
ments listed above, we can apply our previous analysis 
to the cutting lines below the vertex. In addition, how- 
ever, one can ask how the lines continue above the 2 — > 4 
gluon vertex: we show two examples, one of them contain- 
ing a cut (reggeized) gluon. Concentrating on this two- 
gluon state (i.e. we imagine that we have already summed 
over all possible cutting lines below the vertex V), the 
counting rules are quite different: in contrast to the even- 
signature Pomeron, the gluon is a odd-signature reggeon. 
Consequently, the cut gluon is suppressed w.r.t. the uncut 
gluon by one power in as, and this suppression leads to 
the following hierarchy of cutting lines: the cut between 
the gluons belongs to leading order, the cut through one 
of the two reggeized gluons is suppressed by one power in 
Os, the cut through both reggeized gluons is double sup- 
pressed (order a^). A closer analysis of this question is 
under investigation [222]. 



6.1.4 Conclusions 

Corrections due to multiple interactions seem to be im- 
portant in DIS at small x and low Q^; they are expected 
to play a significant role also in multijet production in 
pp scattering. The study of the AGK rules to pQCD pro- 
vides help in understanding the systematics of multiple 
gluon chains. Results described in this review represent 
the beginning of a systematic analysis. We have listed a 
few questions which require further work. 

As an immediate application, we believe that a quanti- 
tative analysis of multiple scattering at HERA will provide 
a useful input to the modeling of final states at the LHC. 



6.2 Experimental consequences 

Main author H. Kowalski 

Experimentally it is easy to differentiate between 
diffractive and single or multiple inclusive final states since 
diffractive states exhibit large rapidity gaps. The multiple 
inclusive final states should also be distinct from the single 
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Fig. 43. The single gluon-ladder contribution to the total "f*p 
cross section. The blob at the lower end of the diagrams con- 
tains the physics below the scale Qg which separates hard from 
soft physics, whereas the blob at the upper end contains hard 
physics that can be described by pQCD. The dashed line de- 
notes the cut. 



inclusive ones since, at least naively, we would expect that 
in the multiple case the particle multiphcity should be con- 
siderably higher. At low x, however, the relation between 
the number of virtual states excited in the interaction (as 
measured by F2) and the final particle multiplicity cannot 
be straightforward since the growth of F2 with decreasing 
X is faster than the multiplicity increase. This may indi- 
cate that the hadronization mechanism may be different 
from the string picture commonly used in the hadroniza- 
tion procedure of single chain parton showers. The in- 
fluence of multiple scattering on the particle multiplicity 
of the final states should also be damped by the energy 
conservation. The cut through several Pomerons leads 
clearly to more gluons produced in the final state, but the 
available energy to produce particles in the hadronization 
phase remains the same. A detailed Monte Carlo program 
is therefore necessary to evaluate this effect. 

The number of diagrams contributing to the reac- 
tion amplitude increases very quickly with the number of 
Pomerons. For the 3-Pomeron amplitude the gluons can be 
paired in 15 possible ways, shown in Fig. l44l with the exam- 
ples of 0-Pomeron, 1-Pomeron, 2-Pomeron and 3-Pomeron 
cuts. For ?n-Pomerons the number of possible gluon pairs 
and also diagrams is: 

(2m-l)(2m-3)(2m-5)....= (2m - 1)!/(2'"-1(to - 1)!). 

Assuming that all the diagrams for a given multi- 
Pomeron exchange amplitude contribute in the same way, 
the above analysis suggests that the probability for differ- 
ent cuts to contribute should be given by the combinato- 
rial factors. This is the content of the AGK rules which 
were obtained from the analysis of field theoretical dia- 
grams well before QCD was established [216] and which 
relate the cross-section, ak , for observing a final state with 
fc-cut Pomerons with the amplitudes for exchange of m 
Pomerons, F^"^^: 



Cfc 



00 



iJ 



^^ ' kl(m-ky. 



pi™) _ 



(80) 



The same result is also obtained from a detailed analy- 
sis of the Feynman diagram contributions in QCD above 
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Fig. 44. 3-Pomeron contributions to the elastic 7*p amplitude. 
All 15 possible diagrams are shown with some examples of 
Pomeron cuts. 



with the oversimplified assumption that only the symmet- 
ric part of the two-gluon couplings contributes [218]. 



6.3 Multiple Interactions in the Dipole Model 

Main author H. Kowalski 

The properties of the multi-Pomeron amplitude and 
of the cut Pomeron cross-sections can be quantitatively 
studied in a dipole model. Along the lines which were dis- 
cussed in section \^ the 7*p interaction proceeds in three 
stages: first the incoming virtual photon fiuctuates into a 
quark-antiquark pair, then the qq pair elastically scatters 
on the proton, and finally the qq pair recombines to form 
a virtual photon. The total cross-section for 7*p scatter- 
ing, or equivalently F2 , is obtained by averaging the dipole 
cross-sections with the photon wave functions, ipir, z), and 
integrating over the impact parameter, b: 



F2 = 



Q' 



Att^c 



Sr 






^ cPb- 



(81) 



Here ?/;*-0 denotes the probability for a virtual photon to 
fiuctuatc into a qq pair, summed over all fiavors and he- 
licity states. The dipole cross-section is assumed to be a 
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function of the opacity ]7: 
dcTa, 



d^^-n^-^^p^-f: 



(82) 



At small- a; the opacity J7 can be directly related to the 
gluon density, xg{x,iJ?), and the transverse profile of the 
proton, T{b): 



Q^—r^as{^l^)xg{x,^?)T{h). (83) 

JVC 

The parameters of the gluon density are determined from 
the fit to the total inclusive DIS cross-section [224]. 
The transverse profile was determined from the exclusive 
diffractive J/W cross-sections [224]. The opacity function 
fi determined in this way has predictive properties; it 
allows to describe other measured reactions, e.g. charm 
structure function or elastic diffractive J/tf' production. 

For a small value of fi the dipole cross-section, eg. (|82|) . 
is equal to fi and therefore proportional to the gluon 
density. This allows to identify the opacity with the sin- 
gle Pomeron exchange amplitude of Fig. 1431 The multi- 
Pomeron amplitude is determined from the expansion: 
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since the dipole cross-section can be expressed as a sum 
of multi-Pomeron amplitudes [225] in the following way: 



du, 
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(86) 
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The cross-section for k cut Pomerons is then obtained 
from the AGK rules, eq. H8U|) . and from the multi-Pomcron 
amplitude, eq. l|5S|l . as: 
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(87) 
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which leads to a simple expression: 

dak fi'' , ^. 



(88) 



r = 0.1 fm 




b (fm) 



Fig. 45. Examples of h dependence of various cut dipole and 
diffractive cross-sections. 



The diffractive cross-section is given by the difference be- 
tween the total and the sum over all cut cross-sections: 






d^h 
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The cut cross-sections determined in the dipole model 
analysis of HERA data have several interesting proper- 
ties shown in Fig. 03J for small dipoles (r = 0.1 fm) the 
opacity fi is also small, so the single cut cross-section, 
CTi, dominates. This leads to particle production emerging 
only from the one-cut pomeron, which should correspond, 
in the context of e.g. the LUND model, to a fragmenta- 
tion of only one string. For larger dipoles (r — 0.6 fm) 
the dipole cross-section starts to be damped in the middle 
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Fig. 46. F2 and the contributions of k-cut Pomeron processes, Fig. 47. Fractions of single (k=l), multiple interaction (MI) 
F^. and diffraction (D) in DIS. 



of the proton (at 6 « 0) by saturation effects. Therefore, 
the single cut cross-section is suppressed in the middle 
while the multiple cut cross-sections, 0-2, (T3, etc, become 
substantial and increasingly concentrated in the proton 
center. These, fairly straightforward properties of dipoles 
indicate that in the central scattering events the multiple 
scattering probability will be enhanced, which may lead 
at the LHC to substantial effects in a surrounding event 
multiplicity. 

The contribution to F2 from the fc-cut Pomeron ex- 
changes are computed in the analogous way to F2 : 



^2 - 



Att^c 



(fr 
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These contributions are shown, together with F2, as a 
function of x for two representative Q^ values in Fig. 1461 
One finds that multiple interaction contributions, i.e. k > 
2, in the perturbative region, at Q^ = 4 GeV^, are sub- 



stantial. In the typical HERA range of x w 10 



10- 



the k = 2 contribution is around 10% of F2 and the con- 
tributions of higher cuts are also non-negligible. For ex- 
ample, the contribution of the 5-cut Pomeron exchanges is 
still around 0.5%, which means that at HERA, many thou- 
sand events may come from this type of process. FigureEl 
shows the fraction of the multiple interaction processes, 



pMI 



pU=2 



+ F, 



k=3 
2 



.i^fc=4. 



F2 in F2, at the same Q 



values. At Q^ = 4 GeV^ the fraction of multiple scattering 
events is around 14% and at Q^ = 40 GeV^ around 6%, in 
the HERA x region, which indicates that the decrease of 
multiple scattering with increasing Q^ is only logarithmic. 
The fraction of diffractive processes, shown for compari- 
son, is of the same order, and drops also logarithmically 
with Q^. The logarithmic drop of the diffractive contri- 
bution expected in the dipole model is confirmed by the 
data [226]. 

The dipole model provides a straightforward extrap- 
olation to the region of low Q^, which is partly pertur- 
bative and partly non-perturbative. Figure 1481 shows the 
contribution to F2 of fc-cut Pomeron processes and the 
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fractions of multiple interactions and diffractive processes 
at Q2 ^ 0.4 GeV^. 

Note also that, as a byproduct of this investiga- 
tion, the ratio of diffractive and inclusive cross-sections, 
F2 /F2 is found to be almost independent of x, in agree- 
ment with the data and also other dipole model predic- 
tions [153,213,226]. The absolute amount of diffractive ef- 
fects is underestimated, since the evaluation of diffraction 
through AGK rules is oversimplified. It is well known [153], 
that a proper evaluation of diffraction should also take 
into account the qqg contribution which is missing in the 
simple AGK schema. 

Hence, we find that the impact parameter dependent 
dipole saturation model [224] reproduces well the main 
properties of the data and leads to the prediction that mul- 
tiple interaction effects at HERA should be of the order of 
diffractive effects, which are known to be substantial. The 
multiple interaction effects should decrease slowly (loga- 
rithmically) with increasing Q^, similarly to the diffractive 
contribution. 



7 Experimental comparisons 

With the luminosity collected at HERA during the past 
years very precise measurements of the proton structure 
function, F2{x,Q'^), have been performed over a large 
range in the fractional proton energy, x, and in the photon 
virtuality, Q^ . The measurements are now limited by sys- 
tematic errors rather than statistical. Parton density func- 
tions have been obtained mainly by fitting the DGLAP 
equations, evolved from an input scale Qq, to the struc- 
ture functions, measured at some scale Q^ . Especially the 
precision data at low Q^ have provided an important in- 
put to various QCD fit analyses. It was recognized early 
that inclusive measurements, like that of structure func- 
tions, are not very sensitive to the new parton dynamics 
expected to appear in the low x region. Instead evidence 
from such dynamics has to be found from investigations 
of hadronic final states in a phase space region where the 
DGLAP governed evolution is suppressed. Thus, a global 
fit, which also includes data from more exclusive processes, 
would further constrain the PDFs. A problem is that mea- 
surements of the hadronic final states suffer from much 
larger uncertainties than the inclusive structure function 
measurements and therefore measurements of many dif- 
ferent complementary processes are desirable. 

Forward jet production in DIS is expected to be sen- 
sitive to new dynamics and early results indeed showed a 
deviation from the predictions of the LO DGLAP model as 
well as of NLO calculations. However, with the inclusion 
of resolved photon contributions, DGLAP provided the 
same level of agreement as the colour diple model (CDM), 
in which the parton emission follows the same scheme as 
in the new dynamics proposed. Only recent studies of fi- 
nal states with a 'forward jet and two additional jets' give 
the first evidence for parton dynamics in which there is 
additional breaking of the fct-ordering compared to that 
predicted by the resolved photon model. 
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Fig. 48. Left: F2 and the contributions of fc-cut Pomeron pro- 
cesses. Right: Fractions of single (k=l), multiple interaction 
(MI) and diffraction (D) in DIS at Q^ = 0.4 GeV^ 



Dijet data may be used to gain better insight into 
the dynamics of the parton evolution and for extract- 
ing updf's. In the low x region boson-gluon fusion pro- 
cesses are dominating and in the LO DGLAP description 
the gluon and the photon collide head on in the hadronic 
center-of-mass system and thus will be produced back-to- 
back. Deviations from this may arise from additional radi- 
ation and if the parton propagator, entering the hard scat- 
tering process, has significant transverse momentum, such 
that the two partons produced in the hard intreraction are 
no longer balanced in transverse momentum. Thus, the 
two jets produced will not be back-to-back in azimuth. 
A measurement of the azimuthal correlation between the 
two jets should be directly sensitive to the predictions of 
models based on different evolution schemes. 

The flavour composition of the final state can also pro- 
vide important information about the evolution and pro- 
duction mechanisms of partons. This has motivated a mea- 
surement of final states with identified strange particles. 
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Although F2 data can be well described by the ex- 
change of a single gluon ladder, it is unlikely that a single 
chain generates large rapidity gaps, which is the signature 
of diffractive processes. The traditional picture of diffrac- 
tive processes is scattering by the virtual photon against 
a pomeron with a partonic structure. Over the past years 
significant progress in the understanding of diffraction has 
been made at HERA, which has led to a modification of 
this description. Data are much better described assuming 
multi-gluon exchange, where a pair of gluons is the mini- 
mum to create a colour singlet state. The niulti-pomeron 
exchange model provides a natural connection between in- 
clusive scattering, diffractive scattering and multiple scat- 
tering given by different cuts through the ladder diagrams 
according to the so called AGK cutting rules, as discussed 
in section |^ Rapidity gaps between high-transverse en- 
ergy jets have been observed at the Tevatron, at a frac- 
tion that is in good agreement with BFKL predictions. 
Also multiple scattering has been studied at the Teva- 
tron, and found to give significant contributions to the 
final state. In ep-collisions at HERA multiple interactions 
can occur in processes where the exchanged photon in- 
teracts via its parton content. Through the possibility to 
control the fraction of the photon momentum, a;^, entering 
into the scattering process, more systematic investigations 
of underlying events may be performed at HERA over a 
wide energy range. 

In general, measurements of final states provide infor- 
mation about the hard scattering process, parton evolu- 
tion, initial and final state radiation and multiple inter- 
actions. Thus, it is important to measure, as accurately 
as possible, the final states in order to test the theoretical 
models. 

In the following the studies of multiple interactions, 
gaps between jets, forward jets and strange particle pro- 
duction will be discussed in more detail. 



7.1 Multiple interactions at the Tevatron and HERA 

Main author J. Turnau 

Since hadrons are composite objects of quarks and glu- 
ons there is a certain probability that collisions between 
hadrons involve more than one parton interaction i.e. we 
have multiple interactions (MI). As a consequence of the 
strong rise of the parton distribution at low x the proba- 
bility to have MI increases with the collision energy and 
the effect at the Tevatron has turned out to be significant. 
At the LHC the contribution from MI will be even larger. 
In electron-proton collisions at HERA MI may occur in 
processes where the exchanged photon is resolved and in- 
teracts via its parton content. The final state of collisions 
with MI will thus contain the products of the primary 
hard collisions, those of additional soft or semihard par- 
ton interactions, contributions from initial and final state 
radiation and from the beam remnants. All products not 
coming from the primary interaction contribute to the so 
called underlying event (UE). 



Effects of MI will influence the total cross section, the 
inclusive jet cross section, the jet multiplicity, the jet pro- 
file, the jet pedestal (the level of transverse energy outside 
the jets), the transverse energy flow and transverse energy 
correlations, the hadron multiplicity, the multiplicity cor- 
relations and may cause large multiplicity fluctuations. 
Experimental data from HERA and the Tevatron have 
been compared to various theoretical models containing a 
description of MI. 



7.1.1 Monte Carlo models for description of multiple 
interactions 

So far multiple interactions are theoretically not well un- 
derstood. The theoretical description is mainly based on 
QCD inspired models, which assume a hard scattering pro- 
cess superimposed on soft or semi-hard interactions. Var- 
ious models differ in how initial and final state radiation 
is taken into account as well as how the hadronization 
process and the beam remnants are treated. 

HERWIG [227,228] assumes that the UE is a soft colli- 
sion between the two beam "clusters" . The parameters of 
this model are tuned to describe experimental results on 
soft hadron-hadron collisions. Also the strength and fre- 
quency parameters of the secondary interactions are sub- 
ject to tuning. There is a possibility to include multipar- 
ton interactions by employing an interface to the JIMMY 
generator [229,230]. To some extent the formalism that is 
used to describe MI in JIMMY is the same as in PYTHIA 
(see below). 

PYTHIA [231] assumes that each interacting beam 
hadron (or resolved photon) leaves behind a beam rem- 
nant, which does not radiate. In contrast to the original 
HERWIG and ISAJET generators PYTHIA uses multi- 
ple parton interactions to enhance the activity of the UE. 
In the simplest version of the PYTHIA multiple inter- 
action model, the transverse momentum cut-off of the 
hard interactions is lowered to p™'° < p™™. The mean 
number of (semi-) hard interactions is given hy < n >— 
o'partonipT''"') / '^nd, where and IS the non-diffractive part 
of the total cross section. The distribution of the number 
of interactions is not uniquely determined. In the simplest 
approach the fluctuations are calculated from a Poisson 
distribution. In the more sophisticated version the number 
of interactions are given by a Poisson distribution for each 
given impact parameter, where the impact parameter de- 
pendence is given by a double-Gaussian overlap function. 
The number of additional interactions is typically of order 
1 -2 . The parton process with the highest transverse mo- 
mentum in the partonic final state can be calculated by 
the quark/gluon 2—^2 matrix element. Additional parton 
interactions in the event are calculated from perturbative 
gluon-gluon scattering processes. 

Simulations of photon-hadron processes have fre- 
quently been performed using the PHOJET generator 
[232]. PHOJET was designed to simulate, in a consistent 
way, all components which contribute to the total photo- 
production cross section. In contrast to PYTHIA, PHO- 
JET incorporates both multiple soft- and (semi-)hard par- 
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ton interactions on the basis of a dual unitarization scheme 
[233]. 

In their initial investigations of UE [234] CDF used 
the ISAJET Monte Carlo [235], which does not include 
multiple scattering a la PYTHIA or HERWIG. Instead 
the beam jets are added assuming that they are identical 
to a minimum bias event at the energy remaining after the 
hard scattering. However, ISAJET did not describe the 
UE data and has not been used in subsequent analyses. 

Generally speaking, the Monte Carlo models which in- 
clude multiple scattering have enough free parameters to 
describe the most important features of data from HERA, 
the Tevatron and of other data found in the JetWeb 
database [236]. A program to tune the model parameters 
is under way. 



7.1.2 Underlying events at the Tevatron 

In the standard analysis of hard scattering events one mea- 
sures jet cross sections and jet properties, which in gen- 
eral are very well described by QCD Monte Carlo models 
and NLO QCD calculations, provided that jet pedestals 
are properly parameterized. The uncertainty in the UE 
contribution to jet events is actually dominating the sys- 
tematic errors for inclusive jet measurements. In order to 
understand the physics of UE, special studies which go 
far beyond a simple parameterization of the energy flow 
outside the jets, are required. 

The CDF collaboration at the Tevatron has performed 
[234, 237] detailed studies of the structure and properties 
of the underlying event in two complementary analyses 
of Run I data at y/ls) = 1800 and ^(s) = 630 GeV. 
The overall event structure was investigated using global 
variables such as charged particle multiplicities and the 
scalar sum of the transverse momenta of charged particles 
as a function of the leading jet momentum. The sensitiv- 
ity to UE is expected to be the highest in phase space 
regions perpendicular to the direction of the leading jet. 
In the first analysis [234] jets were defined by applying the 
simple cone algorithm to charged particles only. Since the 
lower limit of the jet transverse momenta (scalar px sum) 
was chosen as low as 0.5 GeV, UE could be studied in the 
transition region from minimum bias events to events with 
high transverse momentum jets. In a later analysis [237] 
jets were defined using the cone algorithm on calorimetric 
objects with Et > 15 -20 GeV. As shown in Fig. Ei (left) 
the direction of the leading jet in each event is used to de- 
fine different regions in rj — cj) space : "toward" , "away" and 
"transverse" . The "transverse" region is particularly sen- 
sitive to the UE. In ref. [237] the "transverse" region was 
defined as the area in the rj — <f) plane covered by the two 
cones with radii R — y^{Ariy + {Acj})"^ — 0.7 perpendicu- 
lar to the highest energy jet fFig. 1491 right'). On an event- 
by-event basis the regions of "minimal" and "maximal" 
transverse momentum were defined as the regions contain- 
ing the smallest and largest scalar pT sum of charged par- 
ticles, respectively. Such an investigation of the UE helps 
separating the initial and final state radiation component 
from the "beam remnant" components. It can be argued 



Charged Jet #1 
Direction 




Fig. 49. LEFT: Illustration of correlations in the azimuthal 
angle Z\0 relative to the direction of leading charged jet in the 
event. The regions | A(l> \< 60, | A(t> \> 120 and 60 <| Acp \< 
120 are referred to as "towards", "away" and "transverse". 
Each region covers the same range | Arj \ x | Atf} |= 2 x 120°. 
On an event by event basis the regions "transverse mini- 
mum/maximum" are defined to be the ones containing the 
minimum/maximum transverse momentum. 
RIGHT; The phase space regions, as defined in the analy- 
sis [234], shown in the rj ~ (f} plane, where the "transverse" 
regions are given by cones at ±90° to the leading jet direction. 



that transverse energy in the "minimal transverse" region 
{Prp '™™) is due to multiple scattering while the differ- 



ence in transverse momentum between the "minimal-" and 
"maximal transverse" regions AP^^ — Pj. ^™"-'-^_p^ ''"™ 
is a measure of the hard initial/final state radiation con- 
nected to the primary interaction. The CDF analyses have 
established several basic properties of UE, illustrated in 
Figs ISUlf from [234]) andlFTIffrom [237]) and listed below. 

— In the "transverse" regions most sensitive to UE, the 
average number of charged particles and the average 
charged scalar pT sum grow very rapidly with the mo- 
mentum of the leading jet. At prijet) > 5 GeV an 
approximately constant plateau is observed (Fig. 150(1 . 
The height of this plateau is at least twice that ob- 
served in ordinary soft collisions at the corresponding 
energy. Although models including multiple scatter- 
ings (soft or semi-hard) predict a growth of both the 
average number of charged particles and the average 
charged scalar pT sum at low momenta of the leading 
jet, they are not able to describe the data in this region 
{pTiJet) < 5 GeV). 

— For the leading jet above 50 GeV, Pj, ''"™ is almost 
independent on the momentum of the leading jet which 
is correctly described by HERWIG and PYTHIA. 

— The difference APjP increases slowly. 

— Neither PYTHIA nor HERWIG are able to reproduce 
the Pt distribution of tracks in minimum bias events 
(not shown). 

In summary, the QCD models implemented in the 
PYTHIA and HERWIG Monte Carlo programs are able to 
describe the most important features of the UE from the 
Tevatron data. In both cases the agreement is reached only 
after careful tuning of many parameters, in particular the 
regularization scale of the transverse momentum. Clearly 
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Fig. 50. Data (taken from [234]) on the average number of 
charged particles (pr > 0.5 GeV, | t; |< 1) (top) and the scalar 
Pt sum of charged particles (bottom) in transverse region de- 
fined in Fig. 1491 as a function of transverse momentum of the 
leading charged jet compared with Monte Carlo Models. 



the experimental tests of the predictions from PYTHIA 
and HERWIG concerning correlations and fluctuations in 
the UE will be an important challenge over the coming 
years [238]. 

7.1.3 Underlying event energy at HERA 

At HERA, the interaction of electrons and protons via 
the exchange of a quasi-real photon can result in the pro- 
duction of jets. The photon may interact as a point-like 
particle in so called direct process (Fig. 152b ) or it may in- 
teract via its partonic structure such that a parton carry- 
ing a fraction x-y of the photon momentum interacts with 
a parton in the proton. In resolved processes the photon 
remnant can interact with the proton remnant very much 
like in hadron hadron collisions. The center of mass en- 
ergy in the 7p system extends up to 300 GeV, much be- 
low the reach of the Tevatron. Thus the effects of MI at 
HERA are certainly weaker and more difficult to study. 
On the other hand studies of the photon properties from 
measurements of UE at HERA are interesting and com- 
plementary to the measurements at hadron-hadron col- 
liders. The experimental results presented in this section 
have been published by the HI collaboration [239]. They 
are based on a sample where photoproduction events are 
tagged by detecting the scattered electron and it contains 
3 sub-samples : the minimum bias sample (charged track 
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has been tuned to describe the data. 
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Fig. 52. Examples of LO QCD diagrams for photoproduction 
of inclusive jets in direct (a) and resolved (b) photon interac- 
tions. 



with niomentum> 0.3 GeV-|- reconstructed vertex), the 
high- Et sample (total transverse energy in the pseudo- 
rapidity range -0.8 < 77 < 3.3; Et > 20 GeV) and the jet 
sample ( at least 1 jet with Et > 7 GeV). 

At HERA the amount of energy which is carried by 
the photon remnant can be estimated using the variable 
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where a:^*^** is the fraction of the photon energy carried by 

interacting parton, i?^'^ and Ei!^ are the energies of the 
two jets with the highest transverse energies, and ry-''^*^ and 
^jet2 g^j.g ^jj^gjj. pseudorapidities. The energy of the photon, 
E^, is determined from the energy measured in the elec- 
tron tagger. Fig.ESfa) from ref. [239] shows the transverse 
energy density outside the jets of 2-jet events in the cen- 
tral rapidity region. The data decrease as x^^*'' ^ 1 to 
the level measured in deep inelastic ep scattering events, 
dominated by direct photon processes. The dashed line 
in Fig. 1531 indicates the energy density measured in min- 
imum bias events (for which xi^*'^ is not measurable). At 
small xi^*^^ the energy density increases to the level found 
in hadron-hadron collisions (« 0.3 at the SPS and the 
Tevatron). Both PHOJET and PYTHIA with the MI pa- 
rameters suitably tuned (p™*" depending on choice of the 
the photon pdf) are able to reproduce the data. This type 
of the measurement apparently has no analog in hadron- 
hadron collisions. 

Energy-energy correlations are sensitive measures of 
how the energy is distributed over the available phase 
space and provide important information for the modeling 
of UE. The rapidity correlation /2 is defined as 
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Here Et is the total transverse energy measured in the 
HI calorimeter and the other terms refer to transverse 
energies measured in pseudorapidity bins of size Arj = 
0.22 in the jp cms. The average values were extracted from 
all events in the sample. Fig. [531 (b) shows the rapidity 
correlations from the high Et sample. The data show a 
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Fig. 53. a: The transverse energy density outside jets in the 
central rapidity region, | rj-yp |< 1, of 7p collisions as a func- 
tion of the momentum fraction xi^"^ of the parton entering 
the hard scattering process from the photon side. The data 
(taken from [239]) are compared to models with multiple in- 
teractions (PYTHIA mia, PHOJET) and without (PYTHIA). 
The dashed horizontal line marks the energy density level of 
minimum bias events. 

b: the observed rapidity correlations with respect to the cen- 
tral rapidity of 7p collisions, -q* — 0. The dashed (dotted) his- 
togram represents calculations of the QCD generator PYTHIA 
with (without) interactions of the beam remnants. 



short range correlation around the reference bin 77* = 
and a long range anti-correlation which results from the 
hard scattering process. PYTHIA without MI predicts an 
anti-correlation which is too strong. Adding MI i.e. the 
addition of uncorrelated energy to the event results in a 
correct description of the data. The same conclusion holds 
for an event sample where jets are explicitly required (jet 
sample). 

In summary, the underlying event in photoproduction 
events can be consistently interpreted as the superpo- 
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sition of a hard scattering process and interactions be- 
tween the beam spectators, as modeled by PYTHIA and 
PHOJET. Processes with resolved photons at x-y « are 
found to produce 3.5 times the transverse energy den- 
sity of minimum bias events comparable with that ob- 
served in hadron-hadron collisions at the SPS (UAl) and 
the Tevatron (CDF). Studies of energy-energy correlations 
demonstrate that the additional transverse energy in the 
event is not correlated with the hard scattering process. 
Finally, the contribution of higher order radiation to UE 
can be studied separately using the kinematic quantity x^ 
to switch off the beam remnant interactions. 



7.1.4 Explicit Observations of Double Hard Scattering 

The general signature of multiple parton scattering is an 
increase in the transverse energy flow of the event. How- 
ever, in extreme cases, the transverse energy of a sec- 
ondary interaction is sufficient to produce an additional 
pair of jets. The observation of such events is highly impor- 
tant for several reasons. It is sensitive to the phenomenol- 
ogy of multiple parton interactions and provides direct 
information on the structure of the proton in transverse 
space. It is also important for estimating backgrounds to 
processes producing di-boson (W+VF~, etc.) and boson -|- 
jets at the LHC. 

Double parton scattering (DP) in the simplest model 
produces a final state that mimics a combination of two 
independent scatterings. It is customary [240] to express 
the cross section for this process as a product of the cross 
sections for the individual hard scatterings divided by a 
scaling factor, (Jeff- 



crop 



O'Af^B 



2a, 



eff 



The factor m is unity for indistinguishable scatterings and 
has a value of two when it is possible to distinguish be- 
tween A and B. This formula assumes that the number of 
parton-parton interactions follows a Poisson distribution 
but can also use other distributions e.g. Poisson statistics 
for a given impact parameter [241]. The parameter (Jeff 
describes the spatial distribution of partons [242] e.g. for 
a model that assumes a proton with uniformly distributed 
partons Ce// = 11 mb. 

Events with four or more high transverse momentum 
objects (jets, leptons, prompt photons...) is an obvious 
place to look for signatures of multiple hard parton inter- 
actions, although it should be realized that higher order 
QCD processes, for which no exact QCD calculations are 
available yet, are dominating. Only few searches for dou- 
ble parton collisions at the ISR, the SPS and the Tevatron 
have been performed and the results are not very consis- 
tent [240,243,244]. Recently CDF pubhshed [245] a strong 
signal for double parton scattering. In this analysis a value 
of Ceff = 14.5 ± 1.7_2'3 ^^^ ^^^ extracted from data in 
a model-independent way by comparing the number of 
observed double parton events to the number of events 
with hard scatterings at the separate pp collisions within 



the same beam crossing. This represents a significant im- 
provement over previous measurements and may be used 
to constrain models using a parton spatial density. 



7.1.5 IVIultiple interaction component of the underlying 
event at Tevatron and HERA : summary 

Analyses of hadron-hadron and photon-hadron collisions 
at the Tevatron and HERA have firmly established the 
multiple interaction component of the underlying event. 
Only QCD models which include secondary soft or semi- 
hard scatterings a la [242] (PYTHIA, HERWIG, PHO- 
JET) are able to give a reasonable description of the data. 
The energy Sow of underlying events as measured outside 
leading jets was studied in various phase space regions, ap- 
plying conditions which help to disentangle contributions 
from beam-beam interaction and initial/final state radia- 
tion. At HERA the energy available to the photon beam 
remnant was used as an additional constraint. The gen- 
eral structure of the underlying event is reasonably well 
described by Monte Carlo generators like PYTHIA, HER- 
WIG and PHOJET, but a detailed understanding is still 
missing. Studies of underlying events at HERA are not as 
extensive as those by CDF at the Tevatron and it would 
certainly be of great interest to apply the same analysis 
methods to high energy jp, where x~f provides an addi- 
tional "degree of freedom" . The effects of the transverse 
size of hadronic photon on the underlying event, i.e. the Q^ 
dependence, has not been exploited at all so far. The CDF 
Collaboration has reported a firm observation of double 
hard parton scattering in the 7 -f 3 jets final state and 
has made an estimation of the effective cross section for 
double parton scattering. This fact is of paramount im- 
portance for the phenomenological understanding of the 
underlying event, in constraining the multiple interaction 
models [242,246]. 



7.2 Gaps between jets and BFKL 

Main author G. Ingelman 

The observation [247, 248] at the Tevatron of events 
with a rapidity gap between two high transverse-energy 
(Et) jets provides strong evidence for BFKL dynamics 
in terms of color singlet gluon ladder exchange [249]. As 
illustrated in Fig.[^ the process can be described by elas- 
tic parton-parton scattering via a hard color singlet gluon 
ladder. Since there is no color exchanged, no color fields 
(strings) will be formed in between and hence no hadrons 
produced through hadronization in the intermediate ra- 
pidity region. 

In the high energy limit s/|i| S> 1, where the par- 
ton cms energy is much larger than the momentum trans- 
fer, the amplitude for this diagram is dominated by terms 
'^ [as ln(s/|t|)]" where the smallness of ag is compen- 
sated by the large logarithm. These terms are resummed 
in the BFKL equation, which describes the exchange of 
the whole gluon ladder, including virtual corrections and 
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Fig. 54. Hard color singlet exchange through a BFKL gluon 
ladder giving a rapidity gap between two high-p_L jets. 



reggeization of giuons. When solving the equation nu- 
merically it was found that non-leading corrections are 
very important at the non-asymptotic energy of the Teva- 
tron [249,250]. 

Formulating the results as matrix elements for effective 
2^2 parton scattering processes, they could be imple- 
mented in the Lund Monte Carlo Pythia such that par- 
ton showers and hadronization could be added to generate 
complete events. As shown in Fig. 1551 these reproduce the 
data, both in shape and absolute normalization, which is 
not at all trivial. The non-leading corrections are needed 
since the asymptotic Mueller- Tang result has the wrong 
Et dependence. A free gap survival probability parame- 
ter, which in other models is introduced to get the cor- 
rect overall normalization, is not needed in this approach. 
Amazingly, the correct gap rate results from the complete 
model including parton showers, parton multiple scatter- 
ing and hadronization through Pythia together with the 
soft color interaction model [251,252]. The latter accounts 
for QCD rescatterings [253] that are always present and 
if these are ignored one needs to introduce an ad hoc 15% 
gap survival probability factor. 

Related to this is the new results from ZEUS [254] 
on the production of J/ip at large momentum trans- 
fer t in photoproduction at HERA. The data, shown in 
Fig. 123 agree well with perturbative QCD calculations 
[255], based on the hard scales t and nice, for two-gluon 
BFKL color singlet exchange. As illustrated in Fig.|^ not 
only the simple two-gluon exchange is included, but also 
the full gluon ladder in either leading logarithm approxi- 
mation or with non-leading corrections. Using a running 
as does, however, give a somewhat too steep i-dependence 
compared to the data. The conventional DGLAP approx- 
imation provides a good description in the range \t\ < 



v/v 



where this model [256] is argued to be valid due to 
ordered momenta in the gluon ladder (cf. Fig. 157(1 . How- 
ever, the DGLAP model gives a very weak dependence on 
the energy W, which is in contrast to the observed increase 
of the cross-section with energy as also results from the 
BFKL-based calculations [254]. Altogether this provides 
another evidence for BFKL dynamics. 
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Fig. 55. Fraction of jet events having a rapidity gap in |?7| < 1 
between the jets versus the second-highest jet- St. DO data 
compared to the color singlet exchange mechanism [249] based 
on the BFKL equation with non-leading corrections and with 
the underlying event treated in three ways: simple 3% gap 
survival probability, Pythia's multiple interactions (MI) and 
hadronization requiring a 15% gap survival probability, MI 
plus soft color interactions (SCI) and hadronization with no 
need for an overall renormalization factor. Also shown is the 
Mueller- Tang (MT) asymptotic result with a 11% gap survival 
probability. 



7.3 Jets at small-a; 

Main authors L. Jonsson and A. Knutsson 

In the region of low a;-values the interacting parton 
frequently produces a cascade of emissions before it in- 
teracts with the virtual photon. Due to the strong order- 
ing in virtuality, the emissions of the DGLAP evolution 
are very soft close to the proton direction, whereas BFKL 
emissions can produce large transverse momenta in this 
region. Thus, deviations from the DGLAP parton evolu- 
tion scheme are expected to be most visible in a region 
close to the direction of the proton beam. 

HERA has extended the available region in the Bjorken 
scaling variable, xbj, down to values of xbj — 10^^, for 
values of the four momentum transfer squared, Q^, larger 
than a few GeV^, where perturbative calculations in QCD 
are expected to be valid. 

A measurement of the forward jet production cross 
section at small xbj, as proposed by Mueller and 
Navelet [257-259], has long been regarded as the most 
promising test of perturbative parton dynamics. The idea 
is to select events with a jet close to the proton direction 
having the virtuality of the propagator closest to the pro- 
ton approximately equal to the virtuality of the exchanged 
photon. This will suppress an evolution with strong order- 
ing in virtuality as is the case in the DGLAP evolution. 
The additional requirement that the forward jet takes a 
large fraction of the proton momentum, Xj^t = Ejet/Ep, 
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Fig. 56. Differential cross-section da/d\t\ for the process 7 + 
p -^ J/i! + Y. ZEUS data compared [254] to BFKL model 
calculations using leading log (LL) with fixed a^, and including 
non-leading (non-L) corrections with fixed or running q^ as 
well as with a model based on leading log DGLAP. 




Fig. 57. Diffractive vector meson production at large momen- 
tum transfer as described by perturbative QCD hard color sin- 
glet exchange via two gluons and a gluon ladder in the BFKL 
framework [255]. 



such that Xjet ^ xbj opens up for an evolution where 
the propagators are strongly ordered in the longitudinal 
momentum fraction like in the BFKL scheme. Experimen- 
tally this is realized by demanding the squared transverse 
momentum of the forward jet to be of the same order as Q^ 
and Xjet to be larger than a preselected value which still 
gives reasonable statistics. More exclusive final states, like 
those containing a di-jet system in addition to the forward 
jet (called '2-1- forward jet'), provide an additional handle 
to control the parton dynamics. 



Production of forward jets in DIS The HI experiment 
has measured the forward jet cross section [260] using data 
collected in 1997, comprising an integrated luminosity of 



13.7 pb^^. The proton energy is 820 GeV and the positron 
energy is 27.6 GeV which correspond to a center-of-mass- 
energy of ^/s «300 GcV. 

DIS events are obtained by applying the cuts i?e' > 
10 GeV, 156° < Se < 175°, 0.1 < y < 0.7 and 5 GeV^ < 
Q^ < 85 GeV^, where E'^ is the energy of the scattered 
electron, 9^ the polar angle, and y is the inelasticity of 
the exchanged photon. Jets are defined using the inclu- 
sive fct-jet algorithm [261,262] applied in the Breit-frame. 
A forward jet is defined in the laboratory system as hav- 
ing pt.jet > 3.5 GeV and being in the angular range 
7° < 9 jet < 20°. In order to enhance BFKL evolution it is 
required that Xjet > 0.035 whereas DGLAP evolution was 
suppressed in the single differential cross section measure- 
ment by introducing the requirement 0.5 < Pt /Q^ < 5. 

Another event sample, called the '2+forward jet' sam- 
ple, is selected by requiring that, in addition to the forward 
jet, at least two more jets are found, all of them having 
Ptjet larger than 6 GeV. In this scenario the p^/Q^-cut is 
not applied, due to the limited statistics. 

The forward jet cross sections for single and triple dif- 
ferential cross sections are compared to LO (a^) and NLO 
(aj.) calculations of direct photon interactions as obtained 
from the DISENT program. Comparisons of the inclusive 
forward jet cross sections with the DISENT predictions 
for a di-jet final state are adequate, since the forward jet 
events always contain at least one additional jet due to the 
kinematics. The renormalization scale (/i^) is given by the 
average p^ of the di-jets from the hard scattering process, 
while the factorization scale (/ij) is given by the average 
Pj of all forward jets in the selected sample. 

In the analysis of events with two jets in addition to 
the forward jet, the measured cross sections are compared 
to the predictions of NLOJET-I--I-. This program provides 
perturbative calculations of cross sections for three-jet 
production in DIS at NLO (a^) accuracy. In this case the 
scales fir — fif are set to the average p^ of the three se- 
lected jets in the calculated event. 

The NLO calculations by DISENT [263,264] and NLO- 
3ET++ [265] are performed using the CTEQ6M [266] pa- 
rameterization of the parton distributions in the proton. 



Single Differential Cross Section The measured single 
differential forward jet cross sections on hadron level are 
compared with LO (a^) and NLO (a^) calculations from 
DISENT in Fig. &• In Fig. EHb and c the data are com- 
pared to the various QCD models. 

In Fig. I58h . it can be observed that, at small xbj, the 
NLO di-jet calculations from DISENT are significantly 
larger than the LO contribution. This reflects the fact 
that the contribution from forward jets in the LO sce- 
nario is suppressed by kinematics. For small xbj the NLO 
contribution is an order of magnitude larger than the LO 
contribution. The NLO contribution opens up the phase 
space for forward jets and improves the description of the 
data considerably. However, the NLO di-jet predictions 
are still a factor of 2 below the data at low xbj- The 
somewhat improved agreement at higher xbj can be un- 
derstood from the fact that the range in the longitudi- 
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Fig. 58. The hadron level cross section for forward jet produc- 
tion as a function of xbj as measured by HI [260] compared to 
NLO predictions from DISENT (a) and to QCD Monte Carlo 
models (b and c). The shaded band around the data points 
shows the error from the uncertainties in the energy scales of 
the liquid argon calorimeter and the SpaCal electromagnetic 
calorimeter. The hatched band around the NLO calculations 
illustrates the theoretical uncertainties in the calculations, es- 
timated as described in the text. The dashed line in (a) shows 
the LO contribution. 



nal momentum fraction which is available for higher order 
emissions decreases. 

From Fig. EHb it is seen that the CCFM model (both 
set-1 and set-2) predicts a somewhat harder xbj distribu- 
tion, which results in a comparatively poor description of 
the data. 

Fig. 158b shows that the DGLAP model with direct 
photon interactions alone (RG-DIR) gives results similar 
to the NLO di-jet calculations and falls below the data, 
particularly in the low xbj region. The description of the 
data by the DGLAP model is significantly improved if con- 
tributions from resolved virtual photon interactions are in- 
cluded (RG-DIR-t-RES). However, there is still a discrep- 
ancy in the lowest a;_Bj-bin, where a possible BFKL sig- 
nal would be expected to show up most prominently. The 
CDM model, which gives emissions that are non-ordered 
in transverse momentum, shows a behavior similar to the 
RG DIR+RES model. 



Events with Reconstructed Di-jets in Addition to the 
Forward Jet By requiring the reconstruction of the two 
hardest jets in the event in addition to the forward jet, 



different kinematic regions can be investigated by applying 
cuts on the jet momenta and their rapidity separation. 

In this scenario it is demanded that all jets have trans- 
verse momenta larger than 6 GeV. By applying the same 
Pt,jet cut to all three jets, evolution with strong fct-ordering 
is not favored. The jets are ordered in rapidity according 
to -qfwdjet > r]jet2 > Vjeti > Ve ^ith rje being the rapidity 
of the scattered electron. The cross section is measured 
by HI [260] in two intervals of Arji — r]jet2 ~ Vjeti ■ If the 
di-jet system originates from the quarks qi and 52 (see 
Fig. I59|l , the phase space for evolution in x between the 
di-jet system and the forward jet is increased by requir- 
ing that Arji is small and that A'r]2 = rjfwdjet — Vjet-i is 
large. Arji < 1 favors small invariant masses of the di- 
jet system and thereby small values of Xg (see Fig. 159(1 . 
With Ari2 large, Xg carries only a small fraction of the to- 
tal propagating momentum, leaving the rest for additional 
radiation. 

The directions of the other jets are related to the for- 
ward jet through the Arj requirements. When Ari2 is small, 
it is therefore possible that one or both of the additional 
jets originate from gluon radiation close in rapidity space 
to the forward jet. With Arji large, BFKL-like evolution 
may then occur between the two jets from the di-jet sys- 
tem, or, with both Arji and Ar]2 small, even between the 
di-jet system and the hard scattering vertex. By studying 
the cross section for different Ai] values one can test the- 
ory and models for event topologies where the k± ordering 
is broken at varying locations along the evolution chain. 



nrrmnr 9^ 

p"""''"'^ FORWARD JET 

Fig. 59. A schematic diagram of an event giving a forward 
jet and two additional hard jets. These may stem from the 
quarks (qi and 52) in the hard scattering vertex or gluons in 
the parton ladder. Xg is the longitudinal momentum fraction 
carried by the gluon, connecting to the hard di-jet system (in 
this case gi and §2) • 



In this investigation the same settings of the QCD 
models are used as in sections [7.3l while the NLO three-jet 
cross sections are calculated using NLOJET-I--I-. 

From Fig. (HOI it is observed that NLO three-jet gives 
good agreement with the data if the two additional hard 
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jets are emitted in the central region (zir/2 large). It is 
interesting to note that a fixed order calculation (a^), in- 
cluding the log(l/x)-term to the first order in as, is able to 
describe these data well. However, the more the additional 
hard jets are shifted to the forward region {Ari2 small), the 
less well are the data described by NLO three-jet. A pos- 
sible explanation is that the more forward the additional 
jets go, the higher the probability is that one of them, or 
even both, do not actually originate from quarks but from 
additional radiated gluons. NLO three-jet calculates the 
NLO contribution to final states containing one forward 
jet and two jets from the di-quarks, i.e. it accounts for the 
emission of one gluon in addition to the three jets. Since 
the radiated gluon is predominantly soft it has a small 
probability to produce a jet that fulfills the transverse 
momentum requirement applied in this analysis. This re- 
sults in a depletion of the theoretical cross section in the 
small Ar]2 region, which is more pronounced when Arji 
is also small, i.e. when all three jets are in the forward 
region. Consequently a significant deviation between data 
and NLOJET-I-+ can be observed for such events (see the 
lowest bin in Fig. I^Ub ) . Accounting for still higher orders 
in as might improve the description of the data in this do- 
main, since an increased number of gluon emissions would 
enhance the probability that one of the radiated gluons 
produces a jet which is above the threshold on the trans- 
verse momentum. 

As explained above, evolution with strong /c^-ordering 
is disfavored in this study. Radiation that is non-ordered 
in kj_ may occur at different locations along the evolution 
chain, depending on the values of Arji and Ar]2 ■ In a com- 
parsion to QCD models (these figures are not shown, for 
details see [260]) the following observations where made. 
The colour dipole model gives good agreement in all cases, 
whereas the LO DGLAP models give cross sections that 
are too low except when both Arji and Ari2 are large. 
For this last topology all models and the NLO calculation 
agree with the data, indicating that the available phase 
space is exhausted and that little freedom is left for dy- 
namical variations. 

Furthermore it was seen that the '2 -I- forward jet' sam- 
ple differentiates between the CDM and the DGLAP- 
resolved model, in contrast to the more inclusive samples 
where CDM and RG-DIR+RES give the same predictions. 
The conclusion is that additional breaking of the k± order- 
ing is needed compared to what is included in the resolved 
photon model (see Ref. [260]). 



7.4 Production of neutral strange particles in 
deep-inelastic scattering at HERA 

Main author C. Risler 

In deep-inelastic scattering strange particles can be 
produced either if a strange quark is interacting in the 
hard subprocess, or if strange quark pairs are produced 
during the hadronization process. The production of 
strange particles is sensitive to soft and hard parton radia- 
tion of the initial and final state partons and is thus a com- 
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Fig. 60. The cross section for events with a reconstructed 
high transverse momentum di-jet system and a forward jet 
as a function of the rapidity separation between the forward 
jet and the most forward-going additional jet, Ari2 as mea- 
sured by HI [260]. Results are shown for the full sample and 
for two ranges of the separation between the two additional 
jets, Arjx < 1 and Arji > 1. The data are compared to the 
predictions of a three-jet NLO calculation from NLOJET-f-f. 
The band around the data points illustrates the error due to 
the uncertainties in the calorimetric energy scales. The band 
around the NLO calculations illustrates the theoretical uncer- 
tainties in the calculations. 



plementary approach to small x processes. Other sources 
of strangeness can be the decays of charm hadrons or more 
exotic particles like glueballs, pentaquarks or instantons. 

The production properties of strange particles are not 
yet fully understood nor described by the QCD mod- 
els. Since strange particles are also produced during the 
hadronization process, a measurement of strange parti- 
cle production is also a means to test the universality of 
hadronization in e"'"e~, pp or ep collisions. 

The inclusive production cross sections of strange neu- 
tral particles, namely iiTg-mesons and yl-baryons'^, in deep- 
inelastic ep-scattering at HERA were measured with the 
HI detector [267]. The analyzed data were collected in 
the years 1996 and 1997 at a center of mass energy of 



^ By yl-baryons the A particle and its antiparticle A are re- 
ferred to. 



46 



Jeppe R. Andersen et al.: Small x Phenomenology - summary of the 3rd Lund Small x Workshop in 2004 



9 
8 








+^l^±ttH^^^ 


7 
6 
5 
4 
3 
2 
1 

n 


-i- 


■ DATEN 

MEPS, >.s = 0.2 

CCFM, X.s=0.25 



C5 
la 10 



— 1 



V 

1 track unc. 




. . '. . . 


, . , . 1 . , , , 



10 



0.5 1 1.5 



2.5 3 3.5 
Pt [GeV] 



■ri^ 



^ 



• DATEN 

MEPS, ^5 = 0.2 

. - . CCFM, Xc=0.25 



-0.5 



0.5 




2.5 3 3.5 
ri Pt [GeV] 

Fig. 61. Differential Kg and A production cross sections in 
the laboratory frame. 



300 GeV and with an integrated luminosity of 17.8 pb 
The kinematic region 2GeV^ < Q'^ < 100 GeV^ and 
0.1 < y < 0.6 is investigated, where Q^ is the squared 
momentum transfer and y the inelasticity. This allows for 
probing very low Bjorken-cc, a; > 10 



Kg mesons and 



A baryons are reconstructed via the decay to 7r~7r+ and 
TT~p, respectively. The production of Kg and A is mea- 
sured within the visual range, defined by —1.3 < ry < 1.3 
and 0.5 GeV < pt < 3.5 GeV, where rj is the pseudora- 
pidity and px the transverse momentum in the laboratory 
frame. 

Comparisons of the total Kg and A cross sections with 
models using the Lund string hadronization [268-272] 
show that a lower strangeness suppression factor of As ~ 
0.2 — 0.25 is preferred to the default value of A^ = 0.3. 

The differential cross sections in the laboratory and the 
Breit frame are compared to different model predictions, 
namely the MEPS (matrix element and parton shower) 
model using the RAPGAP event generator [273], CCFM 
[24-26] implemented in the CASCADE program [20,50], 
the color dipole model (CDM) [274-279] using D JANGOH 
[280] and to predictions by the HERWIG [227, 228] event 
generator. The HERWIG prediction for the A cross section 
is normalized to the observed total cross section, since 
HERWIG overestimates the cross section by a factor of 3. 

Fig. E] shows the rj and pt dependence of the Kg 
and A cross sections in the laboratory frame compared 
to the model predictions by MEPS and CCFM, using a 
strangeness suppression factor of As — 0.2 and Ag — 0.25 
in the Lund string model, respectively. The preliminary 
data are shown with statistical and systematic errors; the 
systematic uncertainty of the cross section due to the un- 
certainty of the tracking efhciency is separately shown as 
a grey band. 

The rj spectrum of Kg production can not be reproduced 
by the MEPS model, while CCFM gives a better descrip- 
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Fig. 62. T] dependence of Kg and A cross sections in the 
laboratory frame compared to MEPS and a modified MEPS 
model with only BGF hard subprocess. 



tion. In A production a rise in the forward direction, de- 
fined by the direction of the outgoing proton beam, is 
observed, which is not present in any of the models. The 
Pt distribution of the Kg and the A cross section are too 
soft in MEPS as well as in CCFM. 

The CCFM model yields a better description of the 
rj spectra in the data than the MEPS model. In addition, 
CCFM, in its implementation in the CASCADE program, 
allows only for gluon induced hard subprocesses. 

In Fig.EHthe data are compared to a modified MEPS 
model (BGF), where only boson gluon fusion processes are 
taken into account, while quark induced subprocesses are 
switched off. This modified MEPS model gives a slightly 
better description of the data than the standard MEPS 
predictions. 

The differential Kg and A production cross sections 
are investigated as functions of Xp ~ 2\p\/Q and px in the 
Breit frame. 

The Breit frame can be divided into the target hemi- 
sphere of the fragmenting proton and the current hemi- 
sphere in the direction of the incoming photon, which is 
related to the fragmentation of the current quark. 

In the target hemisphere of the Breit frame (fig. I63|l 
all four models underestimate the Kg and A cross section 
at large Xp and the pr spectra are modeled too softly in 
most of the models. 

Only a small fraction of all Kg and A decays is found in 
the current hemisphere of the Breit frame, leading to large 
statistical errors of the differential cross sections shown in 
Fig. 1^ Within these errors CDM gives the best descrip- 
tion of the Xp and px dependence of the differential Kg 
and A cross sections. 

Concluding one can say that none of the models gives 
a satisfactory description of the observed cross sections 
of neutral strange particle production. In particular the 
simulated transverse momentum spectra are too soft. A 
significant increase of yl-baryon production was observed 
in the region < ?/ < 1.3 in the laboratory frame, which 
is not reproduced by any of the models. 

The comparison with QCD models using the Lund 
string hadronization reveals that a lower strangeness sup- 
pression factor than the standard LEP-value seems to be 
preferred at HERA. 
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and pt dependence of Kg and A cross sections in the target 
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The cross sections and comparisons shown here are 
results of a PhD thesis [267]. 



8 Outlook 

Studies of QCD in high energy ep colhsions are interest- 
ing in themselves, as a highly nontrivial theory due to its 
nonlinear nature with a nontrivial vaccuum. It is also im- 
portant in order to fully understand the background in 



attempts to find signals for physics beyond the standard 
model at the LHC and future high energy colliders. 

For the timclike cascades in e^e~-annihilation, exper- 
imental data are reproduced to an extent beyond expec- 
tations, by a perturbative parton cascade (if only the first 
giuon emission is adjusted to matrix elements) followed by 
a model for the subsequent non-perturbative hadroniza- 
tion. To describe the spacelike cascades in ep scattering 
at high energies poses a much more difficult challenge, 
fcx-factorization and leading order BFKL evolution of- 
fer a qualitative frame of reference at small x, but do 
not give a quantitative description of the experimental 
data. Non-leading contributions are large, and the separa- 
tion between perturbative and non-perturbative effects in 
the timelike cascades is not realised in the corresponding 
spacelike processes. 

The non-leading contributions are essential also for the 
behaviour at asymptotic energies. They give asymptoti- 
cally small corrections to the evolution equation, but not 
to its solution. The leading order equation fixes the so- 



lution to the powerlike form 



(with logarithmic 



corrections), but the power A is affected by the non- 
leading terms, which therefore have a very large effect. 
The perturbative-non-perturbative interplay is important 
in two regimes. Firstly, the random walk in Ink'^, charac- 
teristic for the BFKL evolution chain, extends down into 
the soft regime. This problem is further enhanced by a 
running coupling a^. Secondly, the high gluon densities 
at small x imply that unitarity constraints and satura- 
tion become essential. This means that non-perturbative 
effects are important also at larger k±, where the running 
coupling is small. 

Recent progress, described in this report, includes in 
particular: 

— Extending the fc_L -factorization formalism introducing 
two-scale unintegrated and doubly unintegrated PDFs 
and investigation of the importance of the correct kine- 
matics even at lowest order. 

— The solution to the BFKL evolution at NLO, and the 
NLO 7* impact factor. 

— BFKL dynamics in other fields, exemplified by QQ- 
production and away-from-jet energy flow in e+e~- 
annihilation. 

— Studies of unitarity corrections and saturation via the 
Balitsky-Kovchegov equation. 

— Going beyond leading order in the BK equation, where 
in particular energy-momentum conservation has a 
large effect. 

— AGK cutting rules in QCD, multi-pomeron exchange 
and diffraction. 

— Phenomenological applications and comparisons with 
experimental data. Here studies of forward jet and 
heavy quark production are of particular interest. 

Further work is still needed within all these fields. The 
impact parameter dependence and correlations, as well as 
generalisations to eA collisions, need to be studied. This 
is particularly important to get a better understanding of 
high energy proton-proton collisions. To fully understand 
the dynamics of small- a; physics we need in the future also 
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to be able to combine the different routes listed above in a 
unifying formalism, which can simultaneously account for 
the effects of NLO (and NNLO) contributions and uni- 
tarity and saturation effects including multi-pomeron ex- 
change, pomeron loops and diffraction. 

The detailed understanding of small- a; physics is essen- 
tial for the understanding of the underlying event struc- 
ture observed at Tevatron and which is expected to be 
even more significant at the LHC. Small x physics is an 
important issue on its own right and is important also 
for the understanding of the QCD background for any 
searches. Small x physics is very complicated due to the 
large phase space opened but it offers also the possibility 
to understand the transition from a dilute to a dense sys- 
tem in a systematic way and thus contributes much to the 
understanding of complicated processes in general. 
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